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1. The Scientific Problem 

The scientific problem addressed in this research focused on a study of the tropical 
30-60 day oscillation and explanations for this phenomenon. 

1.1 The importance of low frequency oscillations 

The intraseasonal (low-frequency) oscillations with time scales of more than a 
week have received few dynamical explanations in contrast to midlatitude features 
on daily weather charts, whose three dimensional structure is well documented 
by theoretical and observational interpretations. Part of the difficulty is that, in 
general, low-frequency oscillations, such as the 30-60 day oscillations or the Southern 
Oscillation, tend to be defined vaguely during some periods or tend to be obscured 
by the complex interaction of meteorological parameters with different time and 
space scales. These problems are exacerbated by the relative paucity of reliable 
observations (and, consequently, reliable analysis as well) in the tropical latitudes. 

An important aspect of the variability under consideration is in relation to forc- 
ing at climatic changes as characterized by the interannual variability. A possible 
mechanism that may contribute significantly to such variability is the interaction 
between the intraseasonal oscillation and the Southern Oscillation (Lau and Chan, 
1985). At the same time, the intraseasonal oscillations may be a part of a tropical- 
extratropical interaction mechanism. Knutson, et al. (1986) have noted the occur- 
rence of global modes of oscillations with 30-60 day time scales. The 250 mb zonal 
winds in the Southern Hemispheres extratropics during Northern Hemisphere sum- 
mers, were found to be strongly coherent with tropical OLR and zonal winds. This 
research stems from the belief that analysis of such oscillations and investigation 
of mechanisms that cause their development will lead to a better understanding of 
the causes and predictability of climate variability. 
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1.2 Suggested mechanisms for developing the tropical 30- 60 day os- 
cillation 

Theories that have been proposed to explain the intraseasonal oscillation can be 
divided into two categories. First of these include those theories based on exci- 
tation of low-frequency stationary circulations (i.e., nonzonal propagation) in the 
tropical atmosphere by an internal or external source. The second category includes 
those theories based on the eastward propagation of Kelvin and/or Rossby waves 
maintained by convective heating. Neither provides a full explanation to the origin 
and development of the 30-60 day oscillation and the mechanisms that produce its 
period. 

2. Proposed Research and its Objectives 

This research project aimed at studying the linear and nonlinear aspects of the 
30-60 day oscillations. It was originally proposed to: 

(i) develop and apply linear and nonlinear models of atmospheric circulation; the 
nonlinear model being a sigma-coordinate primitive equation model including 
a moisture budget, convective and nonconvective release of latent heat, and 
radiation processes. The effect of spectral horizontal resolution and the simu- 
lation of the diurnal cycles (see (iii) below) on the development of the modeled 
30-60 day oscillation will be addressed; 

(ii) investigate the nonlinear relationships between zonal asymmetries (especially 
wavenumber m=l) and the zonally symmetric mode (wavenumber m=0) in 
producing the energetics necessary for maintaining the oscillation. Such an in- 
vestigation will also contribute to the possible interrelation between the Walker 
Circulation and the slow-moving patterns. 

(iii) We will compare the model output in terms of simulated OLR and reflected 
solar radiation (RSR) to observed patterns of OLR and RSR. The observed 
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fields will be obtained from analysis of ISCCP B1 data (refer to 3.1 below) 
which will provide information not only on a high spatial resolution but also 
on the diurnal cycle which is superimposed on the 30-60 day oscillation. 

(iv) We will evaluate the role of the Hadley circulation for developing the tropical 
low-frequency oscillation, which has been discussed by Anderson and Stevens 
(1987a). It is important that we reexamine the role of the Hadley cell in 
the linear versions of our relatively more complicated model. The phase and 
amplitude of the low-frequency oscillation as calculated in the linear case will 
be analyzed and compared with those calculated using the nonlinear version 
of the 5-level atmospheric model. 

The research was originally proposed for a three year period but only supported 
for two years. Thus some of these objectives were not reached although the research 
has contiuned beyond the expiration of this grant. 

3. Basic approach to research 

The general strategy adopted in this research is to combine analysis of satellite 
data with model calculations. It was initially proposed to use ISCCP data exten- 
sively to diagnose tropical convection and the latent heating associated with this 
convection. As this research progressed, focus was placed on analysis of microwave 
data to identify the moisture content of the tropical atmosphere and its variabil- 
ity. It is proposed that the association of these data with cloud imagery provided 
by ISSCP provides a better relationship between moisture, cloudiness and latent 
heating. Much of the work during year 1 and year 2 focused on both developing 
methods for retrieving water vapor and the subsequent analysis of these data. This 
work has formed the basis of the M. S. thesis of Mr. Daxren Jackson. 

The modeling approach adopted in this study involved: 


(i) Non-linear general circulation model 
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This model was developed by Dr. Hanna (report attached in appendix) and was 
a 5-level primitive equation model. The progress through the second year of this 
research is described in the attached report in Appendix B. The third year of this 
research would have seen the application of this model to the stated objectives of 
the research. 

(ii) Thermodynamic model of the oscillation 

A new and fundamentally different model of the oscillation evolved from the re- 
search conducted under this project. This model proposes that an intrinsic thermo- 
dynamic oscillation in the coupled ocean- atmosphere system exists that has features 
reminiscent of the observed 30-60 day oscillation. This model is described briefly 
below and represents a substantial portion of the Ph. D. thesis of Mr. Qi Hu. 

4. Results 

4.1 Satellite data analysis and retrieval development 

In order to understand the properties of latent heating and factors that govern 
this heating, research was carried out using satellite measurements of microwave 
emission from which vertically integrated water vapor is derived. There were two 
aspects to this research. 

(i) Analysis of SMMR water vapor and a study of its relation to sea surface tem- 
perature (Stephens, 1990). The results of this analysis are given in the paper 
attached in the Appendix A. This climatological study served the purpose 
of providing our understanding of the mean state of the atmosphere and for 
studying perturbations from this mean state. 

(ii) The retrieval of water vapor from SSM/I observations. A new retrieval scheme 
has been developed (Tjemkes et al . , 1991) and analysis of SSM/I water vapor 
data is continuing. A discussion of the retrieval scheme and its performance is 
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given in the attached paper (Appendix A). 

Several months of SSM/I data have been analyzed in conjunction with both OLR 
data and ECMWF wind data in an attempt to develop a better understanding 
between convection and moisture convergence. This analysis forms the bulk of Mr. 
Jackson’s M. Sc. thesis and it is anticipated that this will be completed by the end of 
the summer of 1991. Fig 1 presents an example of water vapor anomalies and OLR 
fields, in the form of Houvmuller diagrams, for the period from September 1, 1987 
to November 30, 1987 centered about the equator. Positive water vapor anomalies 
are hatched. Companion diagrams of zonal wind and OLR are also presented for 
comparison. 

4.2 Thermodynamic model of the oscillation 

A thermodynamic model of the low-frequency oscillation is proposed. The es- 
sential features of this thermodynamic cycle are highlighted in Fig. 2. The essence 
of the cycle lies in the evaporation of moisture from the warm ocean surface lead- 
ing to increases in precipitable water and to a destabilization of the atmosphere 
through heating from below. This leads to an enhanced convective activity and 
to a decrease in water vapor. With more cloud, evaporation is decreased and the 
atmosphere stabilized itself leading to a cloud clearing and to a repeat in the cycle. 
Tests of this idea were conducted by Mr. Hu as part of his ongoing Ph. D. research. 
Mr. Hu began this research under the guidance of Prof. D. E. Stevens as part of 
this research project and continues this research under the direction of Prof. D. 
Randall. 

4.3 The 5-level GCM 

The model was developed and tested during the first two years of this research 
project and was to be run to produce a 5-year simulation for the third and final 
year of this research. A report on this modeling effort is attached. 
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Figure 1: Time-longitude plots of (a) the precipitable water vapor anomalies, (b) the zonal 
wind and (c) the OLR for the period from 1 September 1987 to 30 November 1987. 






Figure 2: Schematic representation of a hypothetical thermodynamic cycle that is proposed 
as a possible explanation of the 30-60 day oscillation in the tropics. 
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On the Relationship between Water Vapor over the Oceans and Sea Surface Temperature 

Graeme L. Stephens 

Department of Atmospheric Science, Colorado State University, Fort Collins ; Colorado 
(Manuscript received 17 July 1989, in final form 30 January 1990) 

ABSTRACT 

Monthly mean precipitable water data obtained from passive microwave radiometry (SMMR) are correlated 
with NMC-blended sea surface temperature data. It is shown that the monthly mean water vapor content of 
the atmosphere above the oceans can generally be prescribed from the sea surface temperature with a standard 
deviation of 0.36 gcm‘ 2 . The form of the relationship between precipitable water and sea surface temperature 
in the range T $ > 15°C also resembles that predicted from simple arguments based on the Clausius-Clapeyron 
relationship. The annual cycle of the mass of SMMR water vapor integrated over the global oceans is shown 
to differ from analyses of fully global water vapor data in both phase and amplitude, and these differences point 
to a significant influence of the continents on water vapor. Regional scale analyses of water vapor demonstrate 
that monthly averaged water vapor data, when contrasted with the bulk sea surface temperature relationship 
developed in this study, reflect various known characteristics of the time mean large-scale circulation over the 
oceans. A water vapor parameter is introduced to highlight the effects of large-scale motion on atmospheric 
water vapor. Based on the magnitude of this parameter, it is shown that the effects of large-scale flow on 
precipitable water vapor are regionally dependent, but for the most part, the influence of circulation is generally 
less than about ±20% of the seasonal mean. 


1. Introduction 

The role of water in the energy budget of the climate 
system is a topic of growing interest within the climate 
research community and has been proposed as the basic 
theme of GEWEX (the Global Energy and Water Cycle 
Experiment (GEWEX) is a component of the World 
Climate Research Program). At present, however, large 
uncertainties exist in the estimation of the various 
components of the global water budget and, conse- 
quently, in certain important components of the global 
atmospheric and oceanic energy budgets. We now rec- 
ognize that the connections between global hydrology 
and the Earth’s energy budget are crucial to the problem 
of climate change. For instance, the predicted global 
warming due to a C0 2 doubling with water vapor feed- 
back is approximately twice the warming predicted 
without feedback (the so-called fixed relative humidity 
assumption of Manabe and Wetherald 1 967 ) . The dis- 
tribution of water vapor, its transport, and divergence 
are also essential ingredients to our understanding of 
the distribution of solid and liquid water in the at- 
mosphere and therefore crucial to the significant and 
perplexing problem of cloud feedback to climate 
change (e.g., Paltridge 1980; Somerville and Remer 
1984;Roeckneretal. 1987; Stephens etal. 1990, among 
others). 


Corresponding author address : Professor Graeme L. Stephens, 
Department of Atmospheric Science, Colorado State University, Fort 
Collins, CO 80523. 


Water vapor is also important to other physical pro- 
cesses that occur in the atmosphere. Water vapor plays 
a decisive role in the transfer of radiation through the 
atmosphere and is important to the transport and re- 
lease of latent heat. The distribution of latent heat re- 
lease is a topic that has received considerable attention 
over the past decade especially with the burgeoning 
interest in the variability of the atmosphere on both 
inter- and intra-seasonal time scales. More specifically, 
several studies appear in the recent literature on the 
topic of intraseasonal variability and most focus on the 
explicit coupling between hydrology and atmospheric 
dynamics. The results of Gill’s ( 1982) study on moist 
dynamics, for example, provide insight into the role of 
moist processes in the physics of low frequency vari- 
ability. Webster (1983) also noted the importance 
played by hydrological processes in a modeling study 
of monsoonal low frequency oscillations. Lau and Peng 
(1988) developed a self-consistent theory to explain 
intraseasonal oscillations of the tropical atmosphere in 
which the effects of latent heating are parameterized 
directly as a function of the convergence of water vapor. 
In both that study and that of Gill (1982), explicit 
relationships between precipitable water vapor and sea 
surface temperature had to be assumed. Other studies 
(e.g., Emanuel 1987; Neelin et al. 1987) have also al- 
luded to the importance of evaporation from the warm 
ocean as a mechanism for sustaining low frequency 
oscillations of the tropical atmosphere. 

These problems underscore the growing need to un- 
derstand, at the very least, the bulk characteristics of 
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the distribution of water vapor over the oceans, and in 
particular, its relationship to sea surface temperature. 
The importance of water vapor to a variety of atmo- 
spheric processes has been recognized for some time 
(e.g., Starr et al. 1969), but the poor coverage of con- 
ventional radiosonde data over the global oceans has 
both hindered our understanding of the distribution 
and transport of water vapor and at the same time 
highlighted the need for satellite-based measurements 
of water vapor. In response to this need, a number of 
different satellite approaches have been proposed over 
the past two decades to measure atmospheric water 
vapor. For instance, there have been several attempts 
to estimate water vapor using passive infrared tech- 
niques. Shen and Smith ( 1973) employed satellite in- 
frared spectrometer radiation measurements ( SIRS-B ) 
taken from the Nimbus-4 satellite to estimate precip- 
itable water. Retrieval of moisture fields from mea- 
surements made with the HIRS instrument on the 
Nimbus-6 and TIROS-N satellites have been made by 
Hillger and Vonder Haar (1981) and Hayden et al. 
(1981). Prabhakara et al. ( 1979) used the high reso- 
lution spectral measurements obtained from the IRIS 
instrument flown on the Nimbus 4 satellite to study 
the distribution of precipitable water over the ocean. 
Perhaps most promising are techniques based on pas- 
sive microwave remote sensing. For example, studies 
such as those of Staelin et al. ( 1 976), Chang and Wilheit 
(1979), Grody et al. (1980), Njoku and Swanson 
(1984), Prabhakara et al. ( 1982, 1985), and Wentz 
(1983) among others have demonstrated the viability 
of microwave sensing of precipitable water over the 
world’s oceans. 

The aim of this study is to examine the relationship 
between monthly averaged precipitable water and sea 
. surface temperature. The basic idea of relating precip- 
itable water to surface measurements is not new. For 
example, Reitan ( 1963), Smith ( 1966), and Viswan- 
adham (1981) attempted to correlate precipitable water 
to surface due point temperature with ranging degrees 
of success. While these studies apply over land, Liu 
and Niiler ( 1984) and Liu ( 1986) considered the cor- 
relation of precipitable water to surface mixing ratio 
over the oceans in an attempt to derive an estimate of 
monthly mean latent heat flux at the ocean surface. In 
the present study, precipitable water is correlated to 
the sea surface temperature. In general, we expect this 
relationship to be complicated by the influence of both 
large-scale atmospheric circulation features on precip- 
itable water as well as local processes that influence the 
relative difference between evaporation from and pre- 
cipitation to the surface. Unlike the previously men- 
tioned studies, the present study aims to provide some 
idea of the relative importance of large-scale processes 
versus local effects on the monthly averaged water va- 
por over the oceans. The data sources used in this study 
are described in section 2, and a simple theoretical re- 
lationship between precipitable water and sea surface 


temperature is developed in section 3. Comparisons 
between the observed relationships and the theoretical 
relationship are presented in section 4, and departures 
from the bulk relationship are shown to identify large- 
scale circulation influences on water vapor. 

2. Data sources 

The water vapor data used in the analyses described 
in this paper are derived from the algorithm developed 
by Prabhakara et al. ( 1982) based on the measured 
microwave brightness temperatures at 2 1 and 1 8 GHz 
obtained by the Scanning Multichannel Microwave 
Radiometer (SMMR) flown on Nimbus-7. An advan- 
tage of their retrieval scheme, and one especially rel- 
evant to the desired goals of the research described here, 
lies in the formulation of the retrieval in terms of the 
difference of brightness temperature between two ad- 
jacent frequencies. By differencing these temperatures, 
the radiometric effects of sea surface emission, and thus 
any significant sea surface temperature biases in the 
data, are therefore minimized. Similarly, the effects of 
cloud water droplets and precipitation on microwave 
emission and thus on the retrieved water vapor are 
largely mitigated by using the difference between two 
adjacent frequencies. In an assessment of their retrieval 
approach, Prabhakara et al. claim that the integrated 
water vapor content could be obtained with an rms 
error of about 0.25 g cm “ 2 based on comparison with 
radiosonde measurements of water vapor. This quoted 
error is similar to that estimated by Alishouse ( 1983 ) 
who used a different retrieval approach on SMMR 
measurements and is also comparable to the estimated 
error in precipitable water obtained from radiosonde 
data. The SMMR data used in this study are monthly 
averages and span a period of approximately five years 1 
(January 1979 to September 1983) including the period 
of the much studied 1982-83 El Nino event. The data 
are grouped into rectangular areas 3° latitude by 5° 
longitude over the oceans from 75°N to 75°S. Figure 
1 shows the 3X5 ocean grid on which the water vapor 
data are distributed. The data-void blackened areas 
cover approximately 55 % of the total global surface 
area. Because of the one day on, one day off cycle and 
occasional missing or bad data, approximately 1 50 ob- 
servations per month were obtained for each 3X5 grid 
box. Furthermore, the sun synchronous orbit of the 
Nimbus 7 provides for only two sets of observations 
each day that data are taken, and the possibility of a 
diurnal bias in the data due to an inadequate sampling 
of the diurnal cycle cannot be ruled out. 

The monthly averaged sea surface temperature 
(SST) data used here are part of the 2.5° global gridded 


' The reported comparison between radiosonde and retrieved water 
vapor was carried out systematically throughout the entire period. 
No appreciable effects of channel drift are evident in these compar- 
isons (Prabhakara, personal communication). 
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Fig. 1. Map portraying the ocean grid used in the analysis of SMMR water vapor data. 
The data void areas ( blackened ) lie over major land masses and islands. 



SST dataset produced operationally by the National 
Meteorological Center (NMC) by “blending” in situ 
and satellite observations. The in situ data consist of 
ship and buoy observations, and the satellite obser- 
vations are obtained using multichannel SST tech- 
niques based on the AVHRR observations on the 
NOAA polar orbiting satellites (e.g., McClain et al. 
1985). Because the first available month of these data 
is May 1979, a total of 52 months of concurrent SST 
and water vapor data are available for the analyses pre- 
sented in this paper. A problem of NMCs SST analysis 
concerns the best way to blend the satellite and ship 
datasets to maximize the information in each. The 
method adopted at NMC was to use two days of in situ 
and satellite data to form the blended product. These 
products were then averaged to produce the monthly 
mean SST data used in this study. Unfortunately the 
two-day period chosen for the blending left significant 
areas of ocean undersampled, and SST information 
for these areas tended to be weighted towards clima- 
tology. As a result, the dataset for the time period that 
coincides with the SMMR data suffers from a tendency 
to smooth over SST anomalies (Gemmil, private com- 
munication). Recognition of these problems helped 
prompt NMC to undertake a modified analysis of SST 
observations by blending data collected over a period 
longer than 2 days (e.g., Reynolds 1988). Unfortu- 
nately, these improved SST diata are not available for 
the SMMR time frame. However, the impact of these 
limitations on the analyses described below are likely 
to be relatively small given the interest of this study on 
climatological relationships between SST and water 
vapor. 


A detailed discussion of water vapor distribution has 
been provided in earlier studies of Chang et al. ( 1984) 
and Prabhakara et al. ( 1985 ). The seasonal climatolo- 
gies of the SMMR water vapor are well described by 
Prabhakara et al. ( 1985). The seasonal climatologies 
of the SMMR water vapor are well described by Pra- 
bhakara et al. ( 1985 ), and two examples are shown in 
Figs. 2a and 3a for the periods December to February 
(DJF) and June to August (JJA) averaged over the 
three non-El Nino years. The respective climatologies 
of SST that match these global water vapor distribu- 
tions are shown in Figs. 2b and 3b. A principal feature 
to note is the large area of atmospheric moisture greater 
than 5 g cm* 2 over the equatorial Indian and western 
Pacific oceans (shaded area) which roughly corre- 
sponds to the pool of water warmer than about 29 °C 
(referred to as the western Pacific warm pool). This is 



Fk3. 2a. The climatological distribution of SMMR precipitable 
water obtained as an average of the three DJF seasons for the years 
indicated. Contours have units of gem' 2 . 



JOURNAL OF CLIMATE 


Volume 3 



1 


637 



Fig. 2b. The climatological distribution of SST obtained as an 
average of the same DJF seasons used to obtain Fig. 2a. Contours 
have units of °C. The warm pool regions exceeding 28°C are high- 
lighted. 


a persistent feature of both seasonal climatologies al- 
though the position and size of the moist air varies 
with season. Also of note are the regions of relatively 
low water vapor located over the eastern subtropical 
waters and associated with cold upwelling off the coasts 
of California, Chile, west Africa, and Angola. 

In contrast to these climatologies, water vapor and 
SST anomalies during May 1983, defined relative to 
the average of the three May months during the period 
1 979 to 1 98 1 , are shown in Figs. 4a and 4b respectively. 
A distinct feature of the water vapor anomaly distri- 
bution is the large area of positive anomaly (in excess 
of 1 .2 g cm ~ 2 ) in the equatorial east Pacific associated 
with a warm SST anomaly in this region. Another ob- 
vious feature is the band of dry air across the Pacific 
that branches out from the equatorial west Pacific east- 
ward and northward into midlatitudes. A less extensive 
negative anomaly also occurs in the south Pacific off 
the coast of Chile. 

The qualitative comparison of the anomaly fields 
provided by Figs. 4a and 4b reveals that a number of 
broadscale features of the water vapor anomaly are re- 
lated to anomalies in SST. The negative water vapor 
anomalies just noted more or less match the negative 
SST anomalies. However, the correlation between wa- 
ter vapor and SST does not always hold. For instance, 
the area of negative SST anomaly in the equatorial 



Fig. 3b. As in Fig. 2b but for JJA. 


west Pacific is not realized as a negative anomaly in 
water vapor. Furthermore, the area of positive SST 
anomaly in the eastern Pacific is less extensive than 
the positive water vapor anomaly in this region. 

3. Analysis procedure 

Despite the complicating influences of both dynam- 
ics and thermodynamics on water vapor distribution, 
it is generally believed that its vertical distribution as- 
sumes a simple and predictable character. With this in 
mind, it is proposed that the vertical profile of specific 
humidity q has the form used previously for example 
by Smith (1966) and Sellers (1973) among others, 
namely 

Q = Qo(PIPo) x ( I ) 

where p is pressure and q 0 ,p 0 are the respective surface 
values. The interpretation of X is such that the scale 
height of water vapor is H/X where H is the atmo- 
spheric scale height. If we take 7 km for a typical value 
of H and 2 km for the water vapor scale height, then 
a typical value of X is 3.5. 

It follows from the integration of ( 1 ) from p = 0 to 
p = p 0 that the precipitable water is 



Fig. 3a. As Fig. 2a but for the three JJA seasons 
of the yean indicated. 


Fig. 4a. The May 1983 precipitable water anomaly distribution 
(contours have units of g cm ' 2 ). 
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where e Q in this expression is the surface vapor pressure 


= PoQo 
° 0.622 


re’ 


(3) 


and r is the surface value of relative humidity. In (3), 
is the saturation vapor pressure (in mb) which is 
determined from an approximation to the Clausius- 
Clapeyron equation of the form 

e* = 17.044e a<r, “ 288) (4) 

where T s is the SST (Kelvin) and a » 0.064 K _1 . 
From (2), (3), and (4), the relationship between pre- 
cipitable water w (in g cm' 2 ) and SST then follows as 


w = 10.82 


r 

l + X 


e a{T,~ 288) 


(5) 


where the factor r/( 1 + X) is deduced below by least 
squares fitting of SMMR derived values of w and the 
NMC SSTs to (5). The relationship expressed by (5) 
is provided only as a convenient point of reference for 
the analysis of the observations described below, and 
no claim is made that ( 5 ) actually represents the ob- 
served relationship between w and SST although it is 
shown that (5) fits the observations well for SSTs in 
excess of 15°C. 


4. Results 

a. Global relationships 

Figures 5a and 5b provide graphical examples of the 
relationship between the seasonal means of SMMR 
derived precipitable water and SST averaged for the 
five June to August (JJ A) seasons and the four Decem- 
ber to February (DJF) seasons contained in the 52- 
month dataset ( upper panels). The points represented 
by open circles are the average of all observations that 
fail within a l °C temperature range centered about the 
given temperature. The extent of the shading above 
and below each of these average points represents the 
standard deviation for the data of each bin. Shown in 
the lower panel are the number of observations con- 
tained in the bin which were used to determine both 
the average and standard deviations. Apparent from 
both diagrams is the expected but dramatic increase in 




Fig. 5. Graphical representations of the relationship between pre- 
cipitable water and SST ( upper panels) obtained from the composite 
of all (a) DJF and (b) JJA seasons contained in the dataset. The 
shading above and below each point is one standard deviation above 
and below the averaged data point. The number of observations used 
to form this average and standard deviation are shown as histograms, 
in the lower panel. 


atmospheric water vapor as SST is increased above 
about 20°C. This feature is also apparent in Fig. 6 
which shows the annually averaged relationship derived 
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Fig. 6. As in Figs. 5a,b but for the annual mean. 


from the composite of all 52 months of data. The stan- 
dard deviation of precipitable water averaged over the 
range of SSTs indicated on this diagram is ±0.36 g 
cm -2 . The similarity of both seasons to the annually 
averaged relationship is further emphasized in Fig. 7. 
The mean DJF and JJA values are shown on this dia- 
gram together with the annual average relationship 
(solid curve). Also included for comparison are three 



Fig. 7. The annually averaged relationship extracted from Fig. 6 
(heavy solid line) compared to the DJF and JJA values from Figs, 
5a,b (symbols). Three examples of the Clausius-Clapeyron expression 
(5) are also shown for the stated values of r/( 1 + X). 


examples of the relationship expressed by (5) with rf 
(1 + X) = 0.1, 0.2, and 0.3. The relationship between 
w and T S1 especially for 5 < T s < 15°C, is shown to 
differ from the simple expression derived from the 
Clausius-Clapeyron relationship using a single value 
of r/(l -F X). The apparent increase of water vapor 
above that predicted from simple theory over this tem- 
perature range occurs at those latitudes marked by sig- 
nificant moisture convergence through meridional 
transports (e.g., Peixoto et al. 1978). However, for 
temperatures T s > 15°C a value of r/( 1 + X) = 0.162 
and a = 0.0686 in (5) provides a least squares fit to 
the data with a standard error of 4.8%. 

Figure 8 presents time series of the monthly mean 
water vapor integrated over the SMMR ocean grid from 
75°N to 75°S( upper panel), over the Northern Hemi- 
sphere (middle panel), and Southern Hemisphere 
( lower panel). The quantity represented by the dashed 
curve is the total integrated water vapor (in kg) or, 
alternatively, averaged precipitable water (in g cm ~ 2 ), 
and the averaged SST is given by the solid curves. The 
mass of water vapor over the global oceans has a mean 
of 7.2 X 10 15 kg which is 49% of the average mass of 
water vapor reported by Trenberth et al. ( 1987). The 
difference can largely be attributed to the reduced sur- 
face area coverage of the SMMR ocean grid compared 
to the full global data used by Trenberth et al. The 
annually averaged precipitable water derived from the 
SMMR data is 2.58 g cm ~ 2 which compares well with 
the values of 2.53 and 2.57 g cm* 2 reported by Rosen 



Fig . 8. Time series of SST ( solid ) and precipitable water ( dashed ) 
averaged over the ocean SMMR grid shown in Fig. 1 . The upper 
panel applies to the average over the entire ocean grid, the middle 
panel to the NH and the lower panel to the SR The water vapor is 
also expressed in units of total water vapor mass ( kg and righthand- 
most scale). 
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et al. ( 1979) and Peixoto and Oort ( 1983) respectively 
but is less than the 2.86 g cm" 2 quoted by Trenberth 
et al. 2 although those authors admit to the possibility 
of a slight high bias in their data. 

Whereas the mean SMMR precipitable water rea- 
sonably agrees with previously determined amounts, 
the amplitude of the annual cycles derived over the 
oceans of each hemisphere disagrees with the hemi- 
spheric cycles deduced from global data. The amplitude 
of the annual cycle deduced from SMMR data is large 
and approximately 30% of the average value of the 
precipitable water over the respective oceans of both 
hemispheres. These cycles largely cancel to produce an 
annual cycle for water vapor over the global oceans 
that is significantly smaller in amplitude than for either 
hemisphere. These results are in direct contrast to those 
of Trenberth et al. who find an amplitude in both the 
global and hemispheric averaged water vapor that is 
approximately twice that shown in Fig. 8. The maxi- 
mum of global SMMR precipitable water over the 
oceans peaks around April and is dominated by the 
annual cycle of water vapor over the southern oceans 
in contrast to the data of Trenberth et al., who find 
maxima in the middle of the year in phase with the 
solar cycle over the Northern Hemisphere. While it is 
possible to suspect the representativeness of the Tren- 
berth et al. data over the oceans as a cause for these 
differences, the differences in the annual cycle of water 
vapor over the oceans compared to that derived over 
the globe are likely to be real and thus point to a sig- 
nificant influence of the continents on global water 
vapor. 

The sea surface temperatures averaged over the 
SMMR water vapor grid show the World Ocean to be 
warmer in April, at the end of the solar heating period, 
than in October, at the end of the cooling period. This 
result is consistent with the oceanic temperature records 
of Levitus ( 1 982 ) although the amplitude is somewhat 
more accentuated. The results suggest that the seasonal 
variation of water vapor over the world oceans reflects 
the seasonal changes in the global SST. 

The SMMR measurements presented in Fig. 8 show 
a systematic decrease in water vapor over the 52 
months analyzed despite the quality checks of Prabha- 
kara which were performed throughout this period at 
certain fixed locations. It is known that the SMMR 
antenna temperatures from which the 2 1 GHz bright- 
ness temperatures are derived systematically drifted 
during this period (e.g., Milman and Wilheit 1985). 
This instrument drift perhaps explains the low fre- 
quency trend apparent in Fig. 8 which is deduced to 
be about 1% per year for the ocean averaged quantity. 
The drift, however, is considerably smaller than this 
on regional scales and generally much smaller than the 
intra-annual variability (shown in Fig. 10). While the 


2 The value derived from the ECMWF analyses has recently been 
revised to 2.55 g cm' 2 ( J. Christy, personal communication ). 


longer term variations shown in Fig. 8 do not invalidate 
the conclusions drawn from this study, the results 
highlight the limitations associated with the use of 
SMMR measurements for studying either long-term 
variability or for inferring global scale changes of water 
vapor. 

b. Regional relationships 

While it is shown that the monthly mean water vapor 
over the oceans is, on the whole, a well-behaved func- 
tion of SST, important deviations from this general 
relationship occur on the regional scale. A number of 
important phenomena might produce regional-scale 
departures in the gross relationship described above. 
For example, low-level moisture convergence such as 
occurs in the vicinity of the ITCZ will produce regions 
of above average moisture. Also advection of moist air 
into a relatively dry region can also act to enhance the 
precipitable water of that region. On the other hand, 
advection of dry, cold air produces a drier than average 
atmosphere as do regions of prevailing subsidence, such 
as within the trade wind regimes where dry air exists 
above an inversion. Thus significant deviations of pre- 
cipitable water from that expected from the simple SST 
relationship reflect these large-scale processes. 

To examine these excursions further, time series of 
water vapor and SST averaged over seven regions were 
studied. The regions chosen are labeled A through G 
in Fig. 9 for further reference. Figure 10 presents time 
series of monthly mean SMMR water vapor (solid 
curves) and SST (dashed curves) for each region. The 
comparisons reveal that the annual cycle of water vapor 
on a regional scale generally correlates with the annual 
cycle of SST except in the western Pacific ( region B ). 
There is also the suggestion that the annual cycle of 
water vapor in the equatorial and warm water regions 
(A, B, D, and F) leads the annual cycle of SST in those 
regions. Of special note is the significant drying that 
occurs in region B during the first few months of 1983 
and the related moistening of the eastern Pacific ( region 
D) during the same period. These deviations from the 
normal water vapor content occur during the El Nino- 
Southern Oscillation episode and are consistent with 
our present understanding of that phenomenon. The 
deep convection normally located over region B and 
associated with the ascending branch of the Walker 
circulation shifts to the eastern Pacific over the anom- 
alous warm water in that region and is replaced by 
large-scale subsidence (e.g., Webster 1987). This sub- 
sidence produces a warming and drying of the middle 
troposphere above a lower, humid boundary layer 
thereby indicating a large-scale dynamical influence on 
atmospheric water vapor. 

Figure 1 1 compares time series of the observed and 
predicted water vapor for the same seven regions and 
is presented to emphasize the differences between ob- 
served monthly mean water vapor and values predicted 
from the gross relationship between water vapor and 




Fig. 9. The location of the seven regions (referenced by A through G) chosen to examine 
the regional scale relationships between water vapor and SST. 



SST derived from Fig. 6. The solid, heavy curve is the Following the work of Prabhakara et al. ( 1 979 ), the 
time series of monthly mean SMMR observations parameter (w - vv)/w is introduced to highlight the 
shown in Fig. 10, and the thin line is the water vapor effects of large-scale motion on atmospheric water va- 
deduced from interpolation of the data contained in por. In the definition of this parameter, w is the ob- 
Fig. 6. The extent of the shading about this thin line served precipitable water and vv is the value estimated 
represents one standard deviation of water vapor above given the SST and the relationship graphically por- 
and below this predicted amount and is again deter- trayed in Fig. 6. Thus, regions influenced by moisture 
mined by interpolation from Fig. 6. The observed water convergence and/or moist air advection are charac- 
vapor amounts generally fall within this range for the terized by positive values of this parameter, whereas 
seven selected regions, and the predicted annual cycles regions characterized by negative values suggest the in- 
of water vapor reasonably match those observed in re- fluence of subsidence or dry air advection on water 
gions A, F, and G. However, there are notable and vapor. Because water vapor is concentrated in the first 
important deviations that exist between the observed few kilometers above the ocean surface, the distribution 
and predicted water vapor amounts which suggest an of ( w - vv)/ w also indicates, in a gross way, the spatial 
influence of large-scale processes. For example, the ef- configuration of the maritime boundary layer (Pra- 
fects of moisture advection are seen in regions C and bhakara et al. 1979) and the magnitude of this quantity 
D. The systematic difference between the observed and also indicates the quantitative effects of large-scale cir- 
predicted water vapor amounts in the east Pacific region culation on precipitable water. 

(D) is due to the influence of low level advection of Figures 12a and 12b present the global distribution 
moisture from the Caribbean Sea into the equatorial of (w — w)/w derived for the DJF and JJA seasons, 
eastern Pacific. The significant moistening of region C The values of this quantity, expressed in Figs. 12a and 
during the boreal summer season can be attributed to 12b as a percentage, vary from around -60% to +20% 
eastward advection of moist air that originates from with large areas over the oceans ranging from ±20%. 
the moister regions of the western Pacific associated Features that appear on these maps reflect various 
with the warm Kuroshio Current and transported east- characteristics of the seasonal mean large-scale circu- 
ward along the northern flank of the subtropical high, lation over the oceans. For instance, the areas of large- 
The observed water vapor amounts in region B also scale subsidence associated both with regions of cold 
tend to be larger than predicted during the boreal sum- upweiling and of large scale subsidence along the east- 
mer season especially during 1980 although this en- em branches of the subtropical anticyclonic circulation 
hancement seems to undergo a marked interannual over the North and South Pacific and Atlantic oceans 
variability. The moistening of the western Pacific region are indicated by negative values of ( w - w)/ w. Positive 
can be related to the increased moisture convergence values of this factor are found in the equatorial regions 
associated with the Asian monsoon. Other large-scale where moisture convergence associated with the ITCZ 
influences on water vapor distribution are evident in prevails. Some other outstanding features of these maps 
region E in which the enhanced wintertime drying re- are the regions of almost continuous negative values 
lates to the large-scale subsidence associated with the of ( w — w)/ w found during JJA in the Southern Hemi- 
accelerating equatorial flow along the eastern flank of sphere west of the date line between 20° and 30°S. The 
the subtropical anticyclones. This explanation is also location of this region broadly agrees with the vertical 
consistent with the vertical motion analysis of Oort motion fields calculated by Oort (1983) and is asso- 
(1983) who shows this subsidence to be strongest dur- dated with the descending branch of the Hadley cir- 
ing the winter season in this region. culation. The relatively moist air over the North Pacific 
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Month /Year 

Fig. 10. Time series ofSST (dashed) and SMMR precipitable water (solid) 
averaged over the seven regions highlighted on Fig. 9. 


in JJA is also consistent with the large northward flux 
of moisture by transient eddies, which is thought to 
occur in this region. A dominant feature of the DJF 
map is the dry air east of Japan and a similar but weaker 
dry feature east of the North American continent. 
These features are respectively associated with the cli- 
matological strong high pressure system located in the 
higher middle latitudes over the cold Eurasian conti- 
nent and a weaker, less extensive high pressure system 
over North America. The dry air perhaps results from 


a combination of strong descent over the eastern por- 
tions of these pressure systems which extends eastwards 
over the Pacific and Atlantic oceans and the advection 
of dry air from the cold land mass over the warmer 
waters of these oceans. The existence of these dry areas 
is also consistent with the large latent and sensible heat 
releases that are perceived to occur in these regions 
(Budyko 1982). Both the seasonal and regional vari- 
ability of ( w - w)/ w, evident for example off the east 
coasts of the Asian and North American continents. 
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Fig. 11 . Time series of SMMR precipitable water (heavy solid) and the precipitable water predicted from the 
annually averaged relationship presented in Fig. 6 (thin solid line, middle curve) for the seven selected regions. The 
shading within the two outer thin lines represents the extent of a standard deviation above and below the predicted 
water vapor. 


may perhaps explafi**hy other studies that attempt 5. Summary and conclusions 
to correlate precipitable water to surface parameters 

show varying degrees of success (e.g., Viswanadham A bulk relationship between monthly mean atmo- 
1981). spheric precipitable water and sea surface temperature 
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was established using 52 months of SMMR microwave 
water vapor data and monthly averaged NMC blended 
SST analyses. It was shown that a composite relation- 
ship between precipitable water and SST resembled 
that derived using simple thermodynamical arguments 
when the sea surface temperature exceeds about 15°C. 
It was also demonstrated that the annually averaged 
mass of water vapor integrated ever the global oceans 
deduced from the SMMR observations is in reasonable 
agreement with certain previous estimates. However, 
both the phase and amplitude of the annual cycle of 
water vapor derived from SMMR, when integrated ei- 
ther over the global oceans or individually over the 
oceans of the Southern or Northern Hemisphere, de- 
viate significantly from some previous findings for the 
globe as a whole. The study ofTrenberth et al. ( 1987), 
for example, determined the amplitude of both the 
global and hemispheric averaged water vapor to be ap- 
proximately twice that found in the SMMR oceanic 
data. The maximum in the annual cycle of global 
SMMR water vapor over the oceans peaks around April 
and is dominated by the annual cycle over the southern 
oceans in contrast to the data ofTrenberth et aL which 
yield a maximum in the middle of the year in phase 
with the solar cycle over the Northern Hemisphere. 

Analyses of the SMMR water vapor on a regional 
scale show that the local departures of observed pre- 
cipitable water from that predicted from SST could be 
used as a tracer of large-scale circulation. Regions of 
large-scale subsidence associated with the mean me- 
ridional circulation and with the circulation about the 
subtropical anticyclones are clearly evident, as are re- 
gions of moisture convergence associated with the 
ITCZ and transient eddy activity in the midlatitudes. 

The basic conclusions that can be drawn from this 
study are as follows: 

(i) The monthly mean water vapor content in the 
atmosphere above the oceans can generally be pre- 
scribed from the SST with an averaged standard de- 
viation of 0.36 gem -2 . The form of the relationship 
between precipitable water and SST in the range T s 



Fig. 1 2a. The DJF climatological distribution of ( w - w)fw over 
the SMMR ocean grid. All DJF seasons of the dataset were used in 
deriving this climatology. Contours are uniUess but when multiplied 
' by 100 are the percentage deviation of the predicted water vapor 
from that observed. 



FIG. 1 2b. As in Fig. 12b but for JJA averaged. 


> 15°C also resembles that predicted from simple ar- 
guments based on the Clausius-CIapeyron relationship. 
For these temperatures a value of r/{ 1 + X) = 0.162 
and a = 0.0686 in (5) provides a least squares fit to 
the data with a standard error of 4.8%. 

(ii) The annual cycle of the globally integrated mass 
of water vapor from SMMR over the oceans differs 
from other analyses in both phase and amplitude. It is 
possible that these differences are real and perhaps point 
to ocean-land differences in water vapor with a possible 
significant influence of the continents on global water 
vapor. 

(iii) Monthly averaged water vapor data (denoted 
as w), when contrasted with the bulk SST relationship 
developed in this study (denoted as w), were shown 
to reflect various known characteristics of the time 
mean large-scale circulation over the oceans. The pa- 
rameter ( w — w)/w was introduced both to highlight 
the effects of large-scale motion on water vapor and to 
provide some quantitative measure of these effects on 
the seasonal mean water vapor. Based on the magni- 
tude of this parameter, it was shown that the effects of 
large-scale flow on precipitable water vapor over the 
oceans is regionally dependent but, for the most part, 
the influence of circulation is generally less than about 
±20% of the seasonal mean. 

The analyses contained in this paper, both in terms 
of the bulk relationship between atmospheric water 
vapor and SST and the relationship between water va- 
por and the large-scale circulation of the atmosphere, 
will hopefully provide a useful framework for com- 
parisons between observed water vapor and that pre- 
dicted by atmospheric general circulation models. For 
instance, it would be useful in the context of water 
vapor feedback in climate simulations to compare the 
correlation between modeled precipitable water and 
SST with the correlation found in this study. It is further 
hoped that comparisons of this type might provide in- 
sight into the behavior of water vapor in the atmosphere 
and point to improvements in our ability to model this 
behavior. 
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1.0 Introduction 


Since relatively small changes in the mass of water vapor on all scales produce significant 
changes of cloudiness, an adequate description of water vapor and its variability is essential for un- 
derstanding the hydrological cycle of the atmosphere and hence the energy balance of the earth's 
climate system. Water vapor is the predominant greenhouse gas and provides a major positive feed- 
back mechanism that influences a carbon dioxide induced climate warming (Manabe and Wether- 
ald 1967; Ravel and Ramanathan 1989). Despite its fundamental importance, the global 
distribution and seasonal variability of both water vapor and its transport are largely uncertain and 
water vapor prediction by climate models, including the strength of the water vapor feedback are 
untested. 

These uncertainties in our understanding are primarily due to the current undersampling of 
the present day routine observing network which is biased to the large northern hemisphere land 
masses. This bias underscores the need for a reliable satellite-based approach to the measurement 
of the global distribution of water vapor. A number of different satellite methods have been pro- 
posed over the past decades. For instance, there have been several attempts to estimate water vapor 
using passive infrared techniques. Shen and Smith (1973) employed satellite infrared spectrometer 
radiation measurements (SIRS-B) taken from the Nimbus 4 satellite to estimate columnar water 
vapor. Prabhakara et al. (1979) applied the high resolution spectral measurements obtained from 
the IRIS instrument flown on the Nimbus 4 satellite to the study of the distribution of columnar 
water vapor over the ocean. Retrieval of moisture fields from measurements made with the HIRS 
instrument on the Nimbus 6 and TIROS N satellites have also been made by Hillger and Vonder 
Haar (198 1) and Hayden et al. (1981) and comparison of water vapor retrieved from operational 
sounders with more conventional data is described by Gruber and Watkins (1981), Illari (1989) and 
others. 


Perhaps most promising are techniques based on passive microwave remote sensing for 
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measuring the vertically integrated water vapor amount (the columnar water vapor, CWV). For ex- 
ample, studies such as those of Chang and Wilheit (1979), Njoku and Swanson (1983), Prabhakara 
et al. (1982, 1985), and Wentz (1983) using data collected from the Scanning Multichannel Micro- 
wave Radiometer (SMMR) flown on Nimbus 7 have demonstrated the viability of microwave sens- 
ing of columnar water vapor over the world’s oceans. 

The principal aim of this study is to provide a method for the retrieval of CWV based on 
the microwave observations of the Special Sensor Microwave Imager (SSM/I) which was launched 
in July 1987 as part of the Defense Military Satellite Program on a sun synchronous satellite. The 
SSM/I is a passive microwave radiometer which measures the upwelling microwave radiation at 
four different frequencies 19.35, 22.235, 37.0 and 85.5 GHz. Measurements of both horizontal and 
vertical polarization are made at 19.35, 37.0 and 85.5 GHz. At 22.235 GHz measurements are 
made in the vertical polarization state only. Of particular interest to the present study are the obser- 
vations at 19.35 GHz. At this frequency the instrument has a spatial resolution at the earth surface 
of 70 x 45 km (3 dB range) and a rms error in the brightness temperature of approximately 0.4 K 
(at 300 K). A complete description of the instrument is given by Hollinger et al. (1987). The SSM/ 
I is planned as an integral part of the DMSP effort was well as NASA’s Tropical Rainfall Measure- 
ment Mission (TRMM) planned for the 1990’s (Simpson etal. 1988). Thus the development of re- 
trievals using these measurements are likely to play a key role in establishing future satellite 
derived water vapor climatologies. Specific advantages of the SSM/I over other microwave instru- 
ments, notably the SSMR, are the improved design and better calibration of the radiometer. Futh- 
ermore, the higher temporal resolution and the larger swath width of the SSM/I relative to the 
SMMR makes the SSM/I observations more suitable for short term water vapor variability studies. 

In this paper we develop a new water vapor retrieval scheme that is validated against inde- 
pendent radiosonde observations and also compared to existing microwave retrieval methods. One 
rationale for the development of a new approach concerns the desire to compare different types of 
retrieval schemes so as to gain some appreciation of the uncertainty associated with the specific 
retrieval method itself. In this regard, we compare the CWV retrieved by the new scheme to that 
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obtained from Wentz (1983), Prabhakara et al. (1982) and Chang and Wilheit (1979). Another mo- 
tivation for the development of the scheme presented here concerns the limitations of applying 
techniques developed for other data sets, specifically for the SMMR, to SSM/I data. These issues 
are examined in section 2 of the paper in which the basic radiative transfer theory is outlined. It is 
demonstrated theoretically how a retrieval approach using the 22.235 GHz and the 19.35 GHz 

channels of the SSM/I leads to ambiguous retrieval estimated of CWV greater the about 40 kgm' 
2 

. The radiative transfer description is then followed by an assessment of the accuracy of the ap- 
proach provided in section 3. Comparison of the method with rawindsonde observations is given 
in section 4 followed by presentation of water vapor retrieval for a three month period in 1987. 
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2.0 Microwave radiative transfer through a cloud free atmosphere 


a) Detailed radiative transfer model 


In the microwave region between 10-40 GHz, the absorption and emission of electromag- 
netic radiation is dominated by water vapor and molecular oxygen. The absorption and emission 
of microwave radiation by liquid water is also important at these frequencies but is neglected here 
although it is a subject of further study. Water vapor exhibits a weak spectral line centered at 22.235 
GHz and a set of oxygen spectral lines are present at 60 GHz. The upwelling microwave radiation 
at the top of the atmosphere follows from the solution of the transfer equation for electromagnetic 
radiation through a cloud free plane parallel atmosphere at a single frequency which is given by 
(Chandrasekhar 1960) 


1 dl(z) 
K(z)sec0 dz 


S(z)-I(z), 


(i) 


where I represents the specific intensity (Wm‘ 2 Hz* l sr* ! ) and S the source function. The frequency 
dependence is hereafter understood. All other symbols are defined in Appendix A. The source 
function S describes the gain of radiative energy due to interaction with atmospheric water vapor 
and oxygen, whereas the second term in the right hand side of (1) represents the loss of radiative 
energy. Assuming an atmosphere in local thermodynamical equilibrium (LTE), the source function 
is given by the Planck function which can readily be approximated by the Rayleigh- Jeans approx- 
imation for the range of microwave frequencies considered here. The radiation received by the sat- 
ellite is expressed in terms of brightness temperature (T b) and is given by the integral solution of 
(1) written as 


T* 1 * v 


(**) 


_ T h,v 

~ A b 


h . dT(z\ z.) 

(z 0 ) ^ (Zq» J T ^ * 


( 2 ) 
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where x the total transmission function defined as the product of the transmission due to oxygen 
(t ojc ) and water vapor (x h o ), 

z' 

t(z o- z '> = = ”P-/(V p M + K « p o« )dz '- (3) 

*0 

The first term of the rhs of (2) describes the microwave emission for the lower boundary given by 


*rK v 
1 b 


(z 0 ) = T.e?' , + 


s~s 


(I ~«^ v ) jT(z‘) 


dx(z n , z') 


o» 
dz’ 


■dz' 


(4) 


with e s the surface emissivity, a function of near surface windspeed, sea surface temperature (SST) 
and polarization. The second term in (2) describes the contribution from the atmosphere. Given the 
vertical distribution of water vapor and oxygen, and the frequency dependance of the absorption 
coefficient, the brightness temperature at the top of the atmosphere is shown to involve three height 
integrals which we evaluate using the trapezoidal rule. The absorption coefficients of water vapor 
and oxygen are deducted from the Gross line shape (Waters 1976). The emissivity of a specular sea 
surface is given by the Fresnel equations (Jackson 1975) using the dielectric constants of Klein and 
Swift (1977). The changes of the surface emissivity due to surface wind speed are calculated with 
the empirical model proposed by Pandey and Kakar (1982). 

We use this microwave radiative transfer model to study two existing retrieval methods. 
The first one is that of Prabhakara et al. (1982) which exploits the difference between the SMMR 
measurements at 18 and 21 GHz to retrieve CWV. This method has the advantage that the retrieved 
CWV is relatively insensitive to the exact value of the surface emissivity (Appendix B). A disad- 
vantage of this method is that for values of CWV > 50 kgm~ 2 , the brightness temperature difference 
exhibits a saturation with respect to precipitable water. This is shown in Fig. 1 where we present 
the difference in the 21 and 18 GHz brightness temperatures as a function of CWV calculated with 
the detailed microwave radiative transfer model using 78 different observed atmospheric temper- 
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ature and humidity profiles. The results for these calculation are represented by the solid dots. For 
small values of CWV, the absorption of microwaves at 18 as well as at 21 GHz, depends linearly 
on CWV, and the brightness temperature difference between the two frequencies is also linear in 
CWV. For CWV > 50 kgm' 2 the emission at 1 8 GHz remains is a linear function of CWV, whereas 
the emission at 21 GHz is not. As a result, the brightness temperature increases no longer a mono- 
tonically with increasing CWV, and the retrieval method of Prabhakara et al. becomes ambiguous 
for CWV > 50 kgm' 2 . The saturation effect becomes even more pronounced for brightness temper- 
ature differences formed using the SSM/I channels one at the center of the weak water vapor line 
at 22.235 GHz and the second one at an adjacent frequency such as at 19.35 GHz. Calculations of 
this brightness temperature difference for the 78 observed profiles mentioned above demonstrate 
that the retrieval of CWV from these measurements using the brightness temperature difference 
method like that of Prabhakara et al. (1982), provides unambiguous results when CWV< 40 kgm' 
2 (Fig.l). 

The second retrieval method published in the litterature is the work of Chang and Wilheit 
( 1979). Their retrieval of atmospheric parameters from passive microwave observations makes use 
of brightness temperature ratios of the form 


R = 


T 2X235 ~ T » 

Tkss-V 


(5) 


This approach of Chang and Wilheit has the same advantage as the method used by Prabhakara et 
al. (1982) since (5) is also insensitive to the actual value of the surface emission (Appendix B). 
However, it can be shown using the theoretical model that (5) exhibits a similar saturation effect 
as that shown for the Prabhakara et al. scheme. Figure 2 shows R for the SSM/I frequencies as a 
function of CWV and highlights the reduced sensitivity of R when CWV exceeds 40 kgm* 2 . 

The analyse of both these methods shows clearly the limitations of using the observations 
at 22.235 GHz for retrieving the CWV especially when CWV > 40 kgm* 2 . Thus an alternate meth- 
od is needed to retrieve CWV in the moist tropical atmosphere. We propose here a method based 
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on the polarization difference at 19.35 GHz. In general e h * e v and the microwave radiation field 
at the top of the atmosphere is polarized. The absorption of the radiation by atmospheric water va- 
por tends to depolarize the radiation field received by the satellite. If the polarization state of the 
radiation field at the surface is known, the difference between the horizontal and vertical polariza- 
tion state of the radiation field at the top of the atmosphere provides a measure of the absorption in 
the atmosphere and thus related to CWV. Figure 3 shows the difference between the brightness 
temperature in the horizontal and in the vertical polarization state as a function of CWV calculated 
with our microwave radiative transfer model. This relation exhibit no saturation effect at high 
CWV as experienced by the previous two methods. However, this new scheme suffers from its own 
drawbacks such as an increased sensitivity to the surface emission. In order to invert the relation 
shown in Fig. 3 we developed the simplified microwave radiative transfer model that is now de- 
scribed. 


b) Approximate radiative transfer model 


The height integrals in (2) and (4) can be evaluated analytically if we assume a specific sim- 
ple expression for the vertical distributions of temperature and water vapor. Assuming T(z)= T s 
+yz, the integral of (2) becomes 

T b z t =* T,(l - (1 -e s ) t 2 ) +Y( Z+ +(!-«,) z'x) , ( 6 ) 


with 


z+ = i z £: x(<z ’ z t )dz * 

*0 


(7) 
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= Jz’^(z 0 * z ') dz ’ 


(8) 


To obtain an estimate of the magnitude of z + and z, we assume that the transmission due to oxygen 
is constant, and use a Taylor expansion for the water vapor transmission to obtain 


z + = x o Jz'ic hj0 (z')p(z') sec0 


i # Tr(f 

n = 0 z' 


K hj0 (z") p (z") dz" sec 9 


dz\ 


( 9 ) 


We can simplify (9) by assuming that the water vapor absorption coefficient is independent of pres- 
sure and temperature and further that the vertical distribution of water vapor is exponential. Thus 
(9) becomes 

Z + = £ (- 1 ) n a n Jz’exp(— ) Mexp( — )dz" 1 dz', 
n-0 - p V. p / 


= X HXzjfz ’ exp ^) fexp (^) -exp (“11 dz ’’ 

n * 0 L P ^ P ^ P JJ 


( 10 ) 


where 


a n = 


K h,VPo +1 <" c0 >“ +t 

n! ' 


(ID 


Use of the binomial theorem gives 


n 


= i (- i ) v ; s b U z ' ex p ( 


(-(n-k+ l))z‘ 


-)dz' 


n-0 


lc-0 


( 12 ) 


with 


bi = (•) 


(W?) 


Since n-k+1 > 0 the integral in (12) can be evaluated yielding 
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«• n Z n + ^ / 


-(n-k+l)z t 


[1 + (n - k + 1) z t /z ] 


(13) 


which can be simplified since z p /z t » 1 such that b® = 1 » b* for V ( k > 1 ) . As a result of this 
simplification. 


z* = Y z n * 2 ( ) n n 

n V p<) <»*» 2 


(14) 


Given thata 0 » a n for n > 1 which follows from the substitution of a maximum anticipated value 
for the CWV of 80 kgm' 2 and an absorption coefficient at 19.35 GHz of approximately 2.5 xlO -3 
m 2 kg* 1 , then 


z * - '.AVo^' ’ Vo ( 1 _ V> (15) 

where we took the first two Taylor expansion terms to approximate the water vapor transmission. 
A similar expression for z can be derived. At a single frequency the multiplicative rule for trans- 
mission holds, that is 

x(z 0 ,z t ) = T(z 0 ,z')T(z’,z t ) 

such that z t - z 0 = (z t - z') + (z' - z 0 ) with z' is at some intermediate level. Then z‘ can be ex- 
pressed as 

« *. * „ x" 

z* = x(ZQ,z t ) £ a nJz ex P(“) ( Jexp(“)dz" dz’ (16) 

n * 0 Zq P z' P 

with a„ given by (1 1). Using arguments similar to those for z + , it follows that 


z* 


= *I 


.n + 2 


2’ 


(n+1) 


(17) 
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and after some simplifications, 


z ‘- TZ P K h 1 oPoV ec 0 a,TZ P (l_ V ) ‘ 


( 18 ) 


at 19.35 GHz. Thus substitution of (15) and (18) into (6) yields the expression for the brightness 
temperature at the top of the atmosphere at 19.35 GHz as 


xh. v 

A b 


(z t ) = T s [ 1 - ( 1 - e,’ v ) t‘] + yz_ ( 1 — t, 


k h 2 o^ ^ T ox + ^ 1 e s 


h, v 


>t ). 


( 19 ) 


The first term on the rhs of (19) describes the contribution of the lower boundary to the radiation 
field at the top of the atmosphere, whereas the second term describes the atmospheric emission con- 
tribution. Equation (19) shows that for a complete transparent atmosphere (t= 1) the brightness 
temperature observed by the satellite is equal to the surface emission. For a completely opaque at- 
mosphere, x=0 and the brightness temperature is given by T s + yz p » T s , where we assume that the 
water vapor is concentrated in the lower parts of the atmosphere and the temperature there is close 
to the sea surface temperature. This shows that by increasing the CWV of the atmosphere, the 
shielding of the surface by water vapor represented by the first term on the rhs of (19), is far more 
important than the increase of the effective emission height, represented by the second term on the 
rhs of (19). 

From these equations it can be seen that for an isothermal atmosphere (Y=0) or for a weakly 
absorbing atmosphere (t » 1), the second term in the rhs of (6) vanishes and the radiation received 
by the satellite is given by 

T b (z t ) = T,(l - (1 - e s )x 2 ) (20) 

which is the formula used by both Prabhakara et al. (1982) and Chang and Wilheit (1979). How- 
ever, the atmosphere is not isothermal, and especially in tropical regions, the transmission at 


22.235 GHz is not close to unity. 
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c) Physical retrieval scheme 


Let A represent the difference between the horizontally and vertically polarized brightness 
temperature. According to ( 19), A and T are related through a cubic relation 


A 




(T,-Y 2 p )r v 



( 21 ) 


where r h ’ v = 1 - e h ’ v . Equation (21) can be solved using standard techniques (cf. Bronshtein and 
Semendyayev 1979) yielding three different complex solutions for x from which CWV has to be 
derived and the physical solution decided on. Rather than proceed in this fashion we choose to sim- 
plify this cubic equation in the following way. Lett* represent the reduced temperature defined by 

t-V^O-T). 

N>x 

and replace this reduced temperature with a climatological SST (T c ). Thus (21) transforms into an 
explicit relation between CWV and A, namely 


CWV 


-[log 


>/(T c r v (l-r h /r v )) 


/(K h sec 0 ). 


ox 


( 22 ) 


We demonstrate below that the theoretical error associated with the simplification of (21) in this 
way is smaller than approximately 2.0 kgm* 2 . We consider® this acceptable given the benefits of 
simplicity and computational efficiency. 

The ratio r 11 / r v in (22) follows from (19) as 


F = 


tJ-t 



(23) 


with 
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T = 


T. + ^d-t). 


ox 


(24) 


Wilheit (1979) showed that the sensitivity of F with respect to T s (i.e. dF/dT) is relatively weak 
such that T can be readily approximated by T c . The latter is taken from the sea surface temperature 
data of Reynolds (1988) archived at the National Center for Atmospheric Research. These data are 
a blend of in situ ship measurements and satellite observations. According to Wilheit (1979) and 
Wentz (1983), F is related to the near surface windspeed by F = go +giV, with go calculated using 
the Fresnel equations and gi taken from the surface emissivity model of Padey and Kakar (1982). 
This linear relation between F and V provides a way of retrieving the near surface windspeed need- 
ed for calculating the surface emissivity. 


The final two parameters needed in the estimation of CWV are the water vapor absorption 

coefficient ( K h o ) and the transmission function for oxygen (x Qx ) at 19.35 GHz. The value of these 

two parameters were obtained from a regression of the total atmospheric transmission calculated 

with the more detailed radiative transfer model described in the previous section and the CWV for 

the 78 observed rawindsonde profiles (Fig. 4). The values used in this study are 
0.47733 _ /T \ 0.68643 


4.1533xl0“ 3 x I ^ 


'T rt \ 


m 2 kg'* and 1.004x10“* x 


,-2 , 


C J 




m 2 kg' 1 with To equals 300 K 


for k. at 19.35 and 22.235 GHz respectively and 0.9807 for t ox at 19.35 GHz. 


hjO 
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3.0 Error analysis and sensitivity study 

In general, the CWV retrieved by such an algorithm is a function of T s and the near surface 
wind velocity (V). The variation of the CWV with respect to SST (i.e. dCWV / dT s ) can be written 
as 

dCWV _ dA -'dA 

dT s ~ dCWV dT/ (25) 

while dCWV / dV follows from 

dCWV _ f <it; V ‘<it; 

dV = VdCWvJ dT s ' (26) 

The sensitivity of the brightness temperature at the top of the atmosphere with respect to uncer- 
tainties in T s (dT b / dSST), CWV (dT b / dCWV) and with respect to uncertainties in V (dT b / dV) 
were estimated from calculations using the detailed microwave radiative transfer model described 
in section 2. The calculations were made for the three standard atmospheres of McClatchey et al. 
(1972) with V equals 10 ms* 1 . The results included in Table 1 list values of dCWV / dT, between 
-0.469 and -0.014 kgm^K* 1 which corresponds to the results reported by Prabhakara et al. ( 1 983). 
This low sensitivity means that the error in CWV will vary between -2.0 and -0.05 kgm* 2 (an ap- 
proximate 3% effect) when the difference between the reduced temperature and the climatological 
T s is about 4 K. 

The sensitivity of the retrieved CWV with respect to the near surface wind speed is given 
in Table 1 as 0.5 kgm*^s. Thus the retrieved CWV is also insensitive to the actual value of V and 
only a rough estimate of V within 4 ms* 1 is sufficient to calculate CWV within 2.0 kgm* 2 . This es- 
timate is provided from the linear relation between F and V proposed by Wilheit (1979) and Wentz 
(1983) discussed in the previous section. 

To study the error in the retrieved CWV due to uncertainties in the brightness temperature 



measurements, some additional calculations were performed using the 78 rawind profiles men- 
tioned before. For this sample of observations the radiometric sensitivity with respect to CWV 
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RS 


CWV 


d 

dCWV 


(Tj-T£)/(T s r v (l-r h /r v )) 


( 27 ) 


1 

T s r v (l-r h /r v ) 


A(Tj-T^) 

ACWV 


( 28 ) 


is evaluated as a function of the CWV. A threshold value of 0.004 m 2 kg * for RScwv derives from 
ACWV » 1 kgm' 2 , A (T{) - T£) =» 0.4K (cf. section 2) and a typical value of 100 K for T s r v (1 - r* 1 / 
r v ). The importance of this threshold value becomes clear if we look at a RS CWV at a particular 
value of CWV. If this RS CWV exceeds the threshold, the resolution in the retrieved CWV is better 
than 1 kgm" 2 whereas if RScwv is smaller than the threshold, the resolution of the retrieved CWV 
will be lower than 1 kgm' 2 . 

The sensitivity of the Chang and Wilheit retrieval method in this study adapted to SSM/I 
channels is 


RScwv = 5cwv(( T n^ 3 - T .)/( T I» 5 - T .)) 
whereas the sensitivity of the adapted Prabhakara et al. scheme is 

RS C?WV = dCWV ^ ^ T 22-235 “ T 19.35) ! ) * 


We normalize these sensitivities according to 


RScwv “Thres 
ll RS cwvll = Max -Thres ’ 


( 31 ) 


where max and Thres refers to the maximum RS CW v value and the threshold for the particular re- 
trieval scheme under consideration (these values are listed in Table 2). Figure 5 shows the normal- 
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ized sensitivities as a function 0 f CWV for the 78 observed profiles mentioned above. For all three 

me ods. die "ormai.zed R Scwv decreases with inching CWV. This means that all methods are 

most sensitive at low CWV values and least sensitive at high CWV values. For the present tetrievai 

sc erne the decrease sensitivity , or high CWV values is less than for me other two schemes as 

expected. Because the methods proposed by both Chang and Wilhei, and Prabhakaraeru/. use the 

observations from the same pair of channels, the normalized RS . u 

2 ’ normalized RS CWV -s change sign around 40 

legm . For these two schemes, CWV values up to 40 kgm* 2 have uncenainties less than 1 kgn,-‘ 
due to uncertainties in the brightness temperamres and greater man 1 kgur 2 forCWV values in 
excess of 40 kgm' 2 . For a realistic range of value, of CWV, me normalized RS CWV f or me present 
retrieval scheme is always positive and thus has an uncenainty less man 1 kgo,' 2 for „„ certainIies 
Of 0.4 K in the measured brightness temperature. 

The error in me retrieval which we can associate with liquid water contamination is how- 
evermore severe than me other errors mentioned above, in me presen, removal scheme no distinc 
t.on between liquid or precipitable water is made. Model calculations of PrabhaJtara e,al ( 1983) 
show tha, the etror in the CWV due to a cloud with a typical water content of 0.25 kg™' 2 may be 
as large as 5 kgm'2. The retrieval of liquid content, in addition to CWV, is therefore required to 
avoid this larger error in the CWV. This is a topic of future research but is to be expected that liquid 
water contamination will only be a source of error in tire instamaneous retrieved CWV. When av- 
eraged over a long period (e.g a mondi), the error in the CWV associated with this cloud contam- 
ination is likely to be snail except in areas of persistent thick clouds where it will introduce a bias 
in the CWV distribution. 



4.0 Results 


a> Com P ar ' son with rawsindsonde observations 


sonde observations. aenved from rawind- 

.m— r 

•™.' , 

the 78 profiles mention * observations different than 

K uies mentioned above, were used. 

Figure 6 shows the scatter diagram of CWV r ^t. • a «. 
the rawindsonde ascents. In general the correlation between the ^ ^ 

is good, with a rms difference is 6.42 kgm‘ 2 (Table (3)) 
one which describes the systematic difference 

chm , , ana could be removed if the appropriate model parameters were 

• m a linear [egression analysis the value for the systematic rms is 4 56 kam* 2 h i 

for the unsystematic rms it is 4.5 1 kgm' 2 . The value for the rm • ' 

for previous micrewave CWV ^e, ™ ™ ■«*- than reported 

1979)mobab| — 

of 1 d I . Pa " “ ^ rawu >*onde measurements, due io the enoneous imerpreution 
of cloud liquid water as CWV dn» »k. . *P on 

in CWV in the drv J T' ° f ' — ^val medtod to Canges 

• retrieval T “ P ^“ dd ““^‘"^"»deU„gofU«,« rt ac««mission(e.g.SS W 

S ° VW an<i ' em> " Wind S P" d ret "'vtU). Test retrievals win, the cubic re- 
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lation (21) showed that the simplification of that relation produces no significant errors in the CWV 
retrieval. 

The uncertainties in the surface emission are related to errors in the retrieval of the surface 
wind speed and to errors in the surface emission model. As mentioned above, the constants used 
in the near surface windspeed retrieval are based on calculations performed by the surface emis- 
sivity model proposed by Pandey and Kakar (1982) and there is some doubt about the constant used 
in the relationship between T h,v and V. In particular, the value of the slope here is -8.159 lO^sm' 1 
(at 291 K) which deviated from the value of -9.310 10* 3 (at 291 K) reported by Wilheit (1979). A 
comparison between the retrieved wind speed and independent observation is therefore needed in 
order to retune these parameters. However, the retrieved CWV is only weakly dependent on the 
retrieved near surface wind speed and such tuning is likely only to alter this rms error slightly. Such 
comparisons were considered beyond the scope of the present paper. Uncertainties in the surface 
emission model itself are probably more severe. Unfortunately there are not enough accurate ob- 
servations to improve the model. 

From the quantitative measures listed in Table (3), we notice that the present retrieval meth- 
od systematically overestimates the CWV for CWV < 20 kgm* 2 . This overestimation is probably 
due to the exclusive use of the 19.35 GHz which results in a low sensitivity of the present retrieval 
method to changes in the CWV in the dry atmosphere. The low sensitivity could be overcome by 
using observations at 22.235 GHz, where the water vapor absorption coefficient is large The results 
of a similar comparison for the adopted Chang and Wilheit applied to the SSM/I observations at 
19.35 and 22.235 GHz are shown in Fig 7 and Table 3. Apparently this particular scheme correlates 
much better with the rawindsonde observations below 35 kgm* 2 than the present method. Since this 
method is independent of the actual surface emission (cf. Appendix B) the better correlations sug- 
gests that the large scatter shown in Fig. 6 is probably due to uncertainties in the surface emission. 
The bias at the low end of the diagram shown in Fig. 6 is removed, but the Chang and Wilheit meth- 
od is unable to retrieve CWV in the moist atmosphere. This is in agreement with the results pre- 
sented in section 2. No retrievals are performed with the Prabhakara et al. method because of the 
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significant ambiguity in their method adapted to the SSM/I frequencies at large values of the C WV 
(cf. Fig. 1) 


b) Horizontal distribution ofCWVfor September October November 1987 


Daily CWV maps were obtained after gridding the retrieved CWV for each observation 
into l°x 1° latitude - longitude grid boxes. These daily maps were then averaged to obtain the sea- 
son mean distribution and its standard deviation as shown in Fig- 8. As noted by Stephens (1990), 
the seasonal mean distribution of CWV also largely resembles the seasonal mean SST distribution 
(not shown). For instance the tongue of low CWV just south of the equator in the East Pacific 
matches the tongue of relatively low SST in that region and the extension of the high CWV region 
in the Western equatorial Pacific is also visible in the seasonal mean SST distribution. This agrees 
with the previously reported results of Prabhakara etal. (1982) and Chang et al. (1984). The nota- 
ble exception is the region of high CWV in the North-Eastern part of the Pacific. This particular 
region is associated with the anticyclonic air motion around a persistent zone of high pressure off 
the coast of Oregon transferring moisture from low-into the mid-latitudes. 

The storm tracks in the northern Alantic show more clearly in the map of standard deviation 
than in the CWV map itself. (Fig* 8 lower panel). The standard deviation map highlights regions of 
high temporal variability, such as found in the storm tracks in the mid-latitude, and within the in- 
tertropical convection zone (ITCZ) and the South Pacific convection zone (SPCZ). It is to be ex- 
pected that these regions of high temporal variation are also regions of large horizontal water vapor 
transfer means of transient eddies and the variance maps shown in the lower panel of Fig. 8 perhaps 
provide an indication of where these transports are greatest. 

When integrated over the ocean grid adopted in this study (approximately 62.7% of the to- 
tal earth surface), an average of 27.63 kgm* 2 (8.86 10 15 kg) of columnar water vapor is retrieved. 
For the analysed season, most of this precipitable water can be found in the Northern Hemisphere 
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where the SST is somewhat higher than in the Southern Hemisphere. In total 29.71 kgm' 2 (4.39 
10 15 kg) is found over the Northern hemisphere oceans compared to 25.84 kgm" 2 (4.47 10 15 kg) 
over the Southern Hemisphere oceans. These values correspond well with the values of 7.2 10* 5 
kg reported by Stephens (1990) as derived from SMMR observations over a smaller area, but is 
larger than the global mean values of 25.3 kgm' 2 (12.9 10 15 kg) and 25.7 kgm" 2 (13.1 10 15 kg) re- 
ported by Rosen et al. (1979) and by Peixoto and Oort (1983) respectively derived from rawind- 
sonde observations. This difference may be due to lower CWV’s over land areas as suggested in 
Stephens (1990). 

The monthly mean distribution of CWV and SST are highly correlated as already men- 
tioned. To study the relation between the SST and CWV in more detail, the SST’s derived from the 
Reynolds dataset were correlated with the CWV derived from the SSM/I observations for Septem- 
ber 1987 in a manner developed by Stephens (1990). The results of such a correlation are compared 
to that obtained by Stephens (1990) for SMMR observations and are shown in Ftg. 9. Apparently, 
the CWV SST correlation for the present dataset does not differ significantly from the results re- 
ported by Stephens. Both relations show an increasing CWV with increasing SST as expected. For 
low values of SST, the CWV of the SSM/I appears about 10 kgm" 2 higher than the SMMR CWV 
values in agreement with the results found from the comparison with rawindsonde observations. 
For intermediate and high SST values the difference between the two relation is small. Some of 
this difference can probably be attributed to the limited size of the present dataset. While the 
SMMR relation is derived from 52 months of observations, the present dataset covers only three 
months. 


b. Intercomparison with other microwave retrievals 


The horizontal distribution of columnar water vapor derived from SSM/I observations us- 
ing the adopted Chang and Wilheit scheme for the analysed season is presented in Fig. 10. This 
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diagram exhibits the same general features shown in Fig. 8, although the CWV’s retrieved in the 
tropical regionsare significantly lower than those shown in Fig. 8. This underestimation is also 
clear from the SST-CWV correlation as shown in Fig. 9. Integrated, over the global oceans with 
the adapted Chang and Wilheit method retrieves an average of 23.22 kgm' 2 (7.45 10 15 kg) from 
these SSM/I observations, which is significant smaller than the values reported above. 

In the present retrieval scheme, LWC is erroneously interpreted as CWV, hence it is antic- 
ipated that the CWV values derived from the SSM/I observations using the present retrieval 
scheme might be too high in certain locations. This is likely to be true for instantaneous retrievals 
of CWV, but not for the monthly and seasonal mean values, since it is likely that only persistently 
thick clouds contribute significantly to this bias. This conclusion is supported from the comparison 
with the monthly mean CWV derived from the SSM/I observations using a method described by 
Wentz (1983). The Wentz method uses a closed form of the transfer equation, similar to our sim- 
plified microwave radiative transfer model described in section 2, to calculate the microwave radi- 
ation received by the satellite. The CWV of the atmosphere follows from an iterative technique that 
matches the calculated radiation field to the observed field. Errors in the Wentz values of CWV due 
to LWC contamination should be significantly reduced since LWC is also retrieved from the SSM/ 
I observations. Moreover, the CWV for regions with moderated or heavy rain are excluded in the 
Wentz retrieval. 

The results of this method for September 1987 is correlated with the sea surface tempera- 
ture for that period and compared to the SST-CWV correlation derived from SSM/I observations 
for the same period using the present method. The results of this comparison are presented in Fig. 

1 1 . In the tropical regions, the CWV values in the Wentz data are higher than the present estimates. 
This difference directly contrasts what would be expected from significant LWC contamination in 
the present scheme. Whether these differences in performance are due to our scheme or to inaccu- 
racies in the Wentz retrieval algorithm is not clear. It appears that the Wentz method retrieves too 
much CWV in the tropical regions. Integrated over the global oceans, the Wentz data produces 1 .03 
10 16 kg (= 31.76 kgm’ 2 ) of water vapor over the global oceans, whereas the present method re- 
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trieves 9.23 10 15 kg (= 29.38 kgm' 2 ) for the same period.A global mean approximately 25.5 kgm' 
2 is reported by Rosen et al. (1979) and Peixoto and Oort (1983) from rawind sonde observations 
(mainly over land. To obtain the mean CWV value given by the Wentz retrieval method the CWV 
a mean value of 13.87 kgm' 2 is required over land. 
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8.0 Summary and Conclusion 


An algorithm to retrieve precipitable water content (CWV) over the global oceans from 
SSM/I observations at 19.35 GHz has been described. The basic philosophy is to utilize measure- 
ments of the polarization state or the radiation field at the top of the atmosphere given the polariza- 
tion state at the sea surface, a function of near surface wind speed and SST. It is shown that the 
radiation field at 19.35 GHz is only weakly dependent on the actual SST and near surface wind 
speed in correspondence with results reported by Prabhakara et al. (1983). A radiative transfer 
analysis is presented and shows the significant limitations of the 22.235 GHz water vapor channel 
for remote sensing water vapor in the moist atmosphere with CWV’s in excess of 40 kgm' 2 . It is 
because of this limitation that the present method is developed, but it is not ideal as it suffers from 
its own sensitivities summarized below. The limitations of the SSM/I for sensing tropical water va- 
por should be a concern addressed in NASA’s TRMM program. 

To estimate the polarization state at the sea surface only rough estimations for SST (within 4 K) 
and the near surface wind speed (within 4 ms* 1 ) are sufficient to retrieve CWV within 2 kgm* 2 . The 
SST information is taken from Reynolds (1988), which is a monthly mean data set based on a 
blending of in situ and bouy observations and satellite observations. The near surface wind speed 
is found from an algorithm proposed by Wilheit (1979) and Wentz (1983). It is concluded that the 
present retrieval scheme is in good agreement with the independent observations provided by raw- 
indsondes. The rawind soundings used for this comparison study were selected to be within 1° and 
within the same day of the SSM/I observation. The standard deviation between the two datasets 
was found to be 6 kgm* 2 , which is somewhat larger than previous reported results. Since the present 
retrieval scheme does not distinguish between liquid and precipitable water, large errors in the in- 
stantaneous retrieved CWV especially in the tropical regions (up to 5 kgm* ) can be expected due 
to liquid water contamination in the present retrievals scheme. However, when integrated over one 
month, the error associated with LWC contamination is probably small, since only thick persistent 
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c'o^ 5 C °" mbU,e “ «« conc.usion is supported by a corapaiison ^ 

m0nthly mea " ^ ° fCWVderiV ' d observations usi „ g the present method with the 

SaItle map dmV ' d tom ,he “* "» "« same period of time but using a different 

retrieval scheme whtch distinguishes between LWC as described by Wen, a (1983). From this com- 

partson ,t ,s found that in general the two methods yields the same results except in the tropics 
whet* the Wentz method retrieves about 10 kgnf* more CWV than the present. Thi s is contrary to 
wha, would be expected from LWC contamination. a futute study, the authors plan to couple the 
present scheme into a method for the retrieval of cloud liquid water content 
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Appendix A. Variables and constants 


e s h,v Surface emissivity in the horizontal (h) and 

vertical (v) polarization state 


go 

gl 

Constants in the linear relation between the 
surface reflectivity ratio (3.4) and the near 
surface wind speed 

[sm" 1 ^ 

I 

Specific intensity 

[Wm^Hz^sr 1 ] 

cwv 

r h,v 
r s 

Columnar water vapor 

Surface reflectivity in the horizontal (h) and 

vertical (v) polarization state 

[kgm' 2 ] 

S 

Source function 

[Wm^Hz'V] 

T b 

Brightness temperature 

[K] 

T c 

Climatological SST 

[K] 

T 

Sea surface temperature (SST) 

[K] 

T 

Reduced temperature 

[K] 

t 

Enhanced temperature 

[K] 

V 

Near surface wind speed 

[ms* 1 ] 

Z P 

e-Folding depth of water vapor 

[m] 

Z t 

Top of the atmosphere 

[m] 

z 0 

Surface 

[m] 

Y 

Temperature lapse rate 

[Km‘ l ] 



[m 2 kg* 1 ] 


K Mass absorption coefficient 

9 A " 8le between v ' ni «ll and propagation direction 

P Density 

r h : o As p but for water vapor 

Pox As p but for oxygen 

T Total transmission 

T hjo as * but for water vapor 

T ox as t but for oxygen 

Constants 


[kgm' 2 ] 

[kgm' 2 ] 

[kgm' 2 ] 


t ox ( 19.35 GHz) = 0.9806 
Y=-6.5 10' 3 [Km* 1 ] 
z p = 2.00 10 3 [m] 

k seed = 4.4516 10* 3 (300/ T s )°- 4 7733 „ 19 35 GHz 
k sec0 = 1.0129 10* 2 (300/T,) 1 - 3 at 22.235 GHz 
go calculated from the Fresnel equations 

g, = -89391xl0~ 3 + 8.104xl0” 3 (T J -275) - 2.0184 x10~ 6 (T s - 275) 2 [sm' 1 ] 
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Appendix B. Alternative retrieval methods 

Let the brighness temperature at the top of the atmosphere be given by (20), then the bright- 
ness temperature difference between two adjacent frequencies can be written as 


T„(o,> -T„(u 2 ) = T s (l-r,t 2 ) -T s (l - r 2 t 2 ) 

(B.i) 

=T s f(I 2 -I 2 ) 

(B.2) 

assuming that r t =* r 2 = i 


Wentz (1983) showed that for the SMMR frequencies f = 0 ( ^-) 

1 s 



such that T s r is approximately independent of the surface temperature. Furthermore, if this differ- 
ence is evaluated for the brightness temperature in the vertical polarization state then this differ- 
ence is also insensitive to the actual wind speed, because dT v /dw s =» 0. (cf. Table l) 


From (20) it is also clear that 


_ T b O>,)-T, _ (M 2 

W-T, " r 2 trj ’ (tj 


(B.3) 


again assuming that Tj * r 2 = f . From we see that R as defined by eq. (5) is independent of the 
actual state of the surface emission. 



References 


29 


Bronshtein, I.N. and K.A. Semendyayev. 1979: Handbook of mathematics. Van Nostrand Reinhold 
CO. 118-122. 


Chandrasekhar S„ 1960: Radiative transfer. Dover Publications, Inc, I -52. 


Chang A.T.C. and Wilheit T.T., 1979: Remote sensing of atmospheric water vapor, liquid water and 

wtnd speed at the ocean surface by passive microwave techniques from the Nimbus 5 satellite. Ra- 
dio Sci., 15 5, 793 - 802 


Chang H.D., Hwang P.H. Wilheit T.T., Chang A.T.C., Staelin D.H. and Rozenkranz P. W, 1 9 84: 

Monthly distributions of precipitable water from the NIMBUS 7 SMMR data. 1. Geophys. Res 89 
D4, 5328 - 5334. 


Gruber, A. and Watkins C.D., 1981: Preliminary evaluation of moisture from NOAA polar orbiting 
satellite sounding radiances. J. Appl. Meteor, 20, 450 - 466 

Hayden D.M., Smith W.L. and Woolf, 1981: Determination of moisture from NOAA polar orbiting 
satellite sounding system. Amer. Water Res. Ass., 115 - 466. 


Hillger D.W. and Vonder HaarT.H., 1981: Retrieval and use of Ugh resolution moisture andsta 

btiity fields from NIMBUS 6 HIRS radiances in pre - convective situations. Mon. Wea. Rev,, 109, 
1788 - 1806. 



30 


Hollinger J.P., 1971: Passive microwave measurements of sea surface roughness. IEEE Trans. 
Geosci. Electron., GE-9 3, 165 - 169. 


Hollinger J.P., Lo R., Poe G., Savage R., and Pierce J., 1987: Special sensor microwave/imager us- 
er’s guide. Naval research laboratory. 


Illari L., 1989: The quality of satellite precipitable water content data and their impact on analyzed 
moisture fields. Tellus, 41 A, 319 - 337. 


Jackson J.D., 1975: Classical Electrodynamics. John Wiley & Sons. 821 - 822. 

Klein L.A. and Swift C.T., 1977: An improved model for the dielectric constant of sea water at mi- 
crowave frequencies. IEEE J. Oceanic. Eng., OE-2 1, 104 - 111. 

Manabe S. and Wetherald R.T, 1967: Thermal equilibrium of the atmosphere with a given distri- 
bution of relative humidity. J. Atmos. Sci., 24 3, 241 - 259. 

McClatchey R.A., Fenn R.W., Selby J.E.A., Volz F.E., and Garing J.S., 1972: Optical properties of 
the atmosphere. AFCRL-72-0497. 

Njoku E.G. and Swanson L., 1983: Global measurements of sea surface temperature, windspeed 
and atmospheric water content from satellite microwave radiometer. Mon. Wea. Rev., Ill, 1977 - 


1987. 



31 


Pandey P.C. and Kakar R.K., 1982: An empirical microwave emissivity model for a foam covered 
sea. IEEE J. Oceanic. Eng., OE-7 3, 135 - 140. 


Peixoto J.R. and Oort A.H., 1983: The atmospheric branch of the hydrological cycle and climate. 
in Variations on the Global Water Budget, D. Riedel, 5 - 65, Hingham, Mass. 


Prabhakara C., Chang H.D. and Chang A.T.C., 1982: Remote sensing of precipitable water over 
the oceans from Nimbus-7 microwave measurements. J. Appl. Meteor., 21, 59 - 68. 


Prabhakara C., Dalu G., Lo R.C. and Nath N.R., 1979: Remote sensing of seasonal distribution of 
precipitable water vapor over the oceans and the interference of boundary layer structure. Mon. 
Wea. Rev., 107, 1388- 1401. 


Prabhakara G, Short D.A., Wiscombe W., and Fraser R.S., 1986: Rainfall over oceans inferred 
from Nimbus-7 SMMR: Application to 1982-833 El Nino. J. Climate Appl. Meteor., 25, 1464 - 
1474. 


Prabhakara C., Wang I., Chang A.T.C. and Gloersen P., 1983: A statistical examination of Nimbus- 
7 SMMR data and remote sensing of sea surface temperature, liquid water content in the atmo- 
sphere and surface wind speed. J. Climat. Appl. Meteor., 22, 2023 - 2037. 


Ravel A. and Ramanathan V., 1989: Observational determination of the greenhouse effect. Nature, 
342, 758 - 761. 



32 


Reynolds R.W., 1988: A real time global sea surface temperature analysis. J. Climate, 1, 75 - 86. 


Rosen R.T., Salstein D.A. and Peixoto J.P., 1979: Variability in the annual fields of large-scale at- 
mospheric water vapor transport. Mon. Wea. Rev., 107, 26-37. 


Shen W.C. and Smith W.L., 1973: Statistical estimation of precipitable water with SIRS-B water 
vapor radiation measurements. Mon. Wea. Rev., 101, 24 - 32. 


Simpson J., Adler R.F. and North G.R., 1988: A proposed tropical rainfall measuring mission 
(TRMM) satellite. Bull. Amer. Meteor. Soc., 69, 278 - 295. 


Stephens G.L., 1990: On the relationship between vapor over the ocean and sea surface tempera- 
ture. J. Climate, 3, 634 - 645. 


Waters J.W., 1976: Absorption and emission of microwave radiation by atmospheric gases. In: 
Methods of experimental physics, M.L.Meeks ed. Academic Press Section 2.3. 


Wentz F.J., 1983: A model function for ocean microwave brightness temperature. J. Geophys. Res., 
88 C3, 1892 - 1908. 

Wetherald R.T. and Manabe S., 1988: Cloud feedback processes in a general circulation model. J. 
Atmos. Sci., 45, 1397- 1415. 



33 


WUhei, T.T.. 1979: Tke effeao/.Mcn ,He mi cro»ave 
GHz. J. Geophys. Res., 84 C8, 4921 - 4926. 


emission from the ocean’s surface at 37 


Wilmott C.J., 1982 : Some 
Soc., 63, 1309 - 1313. 


comments on the evaluation of model performance. 


Bull. Amer. Meteor. 



34 


Figure and Table Captions 

Fig. 1: Relation between the monochromatic brightness temperature difference at the top of the at- 
mosphere between two adjacent frequencies as a function of columnar water vapor calculated for 
78 different rawindsonde ascents using our detailed microwave radiative transfer model. The solid 
points represent the calculated difference between the brightness temperatures at 21 and 18 GHz 
whereas the open points represents the brightness temperature difference at 22.235 and 19.35 GHz 


Fig. 2: Relation between R defined by (5) in the text and CWVcalculated for 78 different rawind- 
sonde ascents. 


Fig. 3: The calculated difference between the horizontal and vertical polarized brightness temper- 
ature at the top of the atmosphere as a function of CWV 


Fig. 4: Correlation between the transmission function evaluated by a microwave radiative transfer 
model at 19.35 GHz (solid points) and at 22.235 GHz (open points) for 78 different rawindsonde 
ascents as a function of columnar water vapor. 


Fig. 5: Comparison of the normalized sensitivity of three retrieval methods as a funcion of CWV 


Fig. 6: Comparison of columnar water vapor derived from SSM/I observations over the global 
oceans for September, October and November 1987 using the present retrieval method, with co- 
lumnar water vapor derived from quasi-coincident rawindsonde observations 
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Fig. 7: As Fig. 6 but now for the Chang and Wilheit retrieval method asapted for SSM/I observa- 
tions. 


Fig. 8(a): Horizontal distribution of the seasonal mean columnar water vapor derived from SSM/I 
observations over the global oceans for September, October and November 1987 using the present 
retrieval method. 


Fig 8(b): The standard deviation in CWV indicating departures of instantaneous observations from 
the seasonal mean. These deviations are averaged to produce seasonal mean values for SON 87. 


Fig. 9: Correlation between sea surface temperature and columnar water vapor as derived from 
SSM/I using the present retrieval method (open circles), using the adopted Chang and Wilheit 
method (triangles). As shown is the correlation as derived from SMMR observations (closed cir- 
cles) 


Fig. 10: As Fig. 8 but for the adopted Chang and Wilheit retrieval method 


Fig. 11. Correlation between sea surface temperature and columnar water vapor as derived from 
SSM/I observations for September 1987 using the present retrieval method (open circles), and ac- 
cording to the Wentz dataset (closed circles) 


Table 1: Theoretical radiometric metric sensitivities at 19.35 GHz for a standard tropical (TRP), 
midlatitude summer (MLS) and midlatitude winter (MLW) atmosphere. 
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Table 2: Values for the maximum and threshold radiometric sensitivities of the retrieved CWV with 
respect to observational errors in the brightness temperature 


Table 3: Quantitative measures of a linear regression analysis between SSM/I retrieved CWV and 
quasi-coincident rawindsonde observed CWV. 
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Table 1 


Atmospheric profile 

TRP 

MLS 

MLW 

CWV [kgm' 2 ] 

57.18 

40.96 

10.16 

T,[K] 

300.0 

294.0 

273.0 

dT v /dT s [K/K] 

0.311 

0.264 

-0.291 

dT h /dT s [K/K] 

0.135 

0.099 

-0.298 

dT v /dCWV [K/kgtn 2 ] 

0.568 

0.628 

0.658 

dT h /dCWV [K/kgm' 2 ] 

0.943 

1.044 

1.176 

dT v /dWj [K/ms’ 1 ] 

0.248 

0.295 

0.361 

dCWV/dT s [kgm* 2 /K] 

-0.469 

-0.397 

-0.014 

dCWV/dWj [kgm^/ms' 1 ] 

2.160 

2.168 

2.135 
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Table 2 


Retrieval Algorithm 

Maximum Sensitivity 

Threshold sensitivity 


[mV 1 ] 

[mV 1 ] 


Chang and Wilheit (1979) 

-9.0732 10' 3 

0.005 

Prabhakara et al. (1982) 

6.3994 10’ 3 

-0.005 

Present 

-7.9964 10° 

0.002 



Table 3 


Variable 

Raobs 

Present 

ACW 

Number of observations 

1188 



Mean [kgm* 2 ] 

16.93 


17.61 

Standard deviation [kgm* 2 ] 

11.20 


9.36 

Offset [kgm' 2 ] 



5.44 

Slope 



0.670 

Mean absolute error [kgm' 2 ] 



3.824 

Root mean square error [kgm' 2 ] 



5.995 

Systematic rms [kgm* 2 ] 



4.122 

Unsystematic rms [kgm* 2 ] 



4.353 
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1. INTRODUCTION 

Although eastward propagation is usually regarded as 
an essential feature of the low frequency Madden- Julian 
oscillation observed in the tropical atmosphere (e.g., Chang, 
1977, Lau and Peng, 1987, ana Wang and Chen, 1989), many 
observations indicate that there is an important stationary or 
quasi- stationary component of the oscillation. Yasunari 
(1979), for example, investigated the stationary 30-60 day 
variation in upper tropospheric cloudiness in the Asian 
summer monsoon region, fn a case study of the 30-60 day 
oscillation, Hsu et al. (1990) found that a case may not 
complete the global circle while its phase cycle has been 
complete. Some of their cases show a strong stationary 
oscillation of the divergence, outgoing longwave radiation 
(OLR) and other Helds. A recent observational study by 
Weickmann and Khalsa (1990) offers further evidence that the 
Madden-Julian oscillation has an important stationary 
component. 

In this paper, we present evidence that intraseasonal 
oscillations can be produced by local radiative and convective 
processes. This suggests that the observed propagating 
* 4 Madden-Julian wave’ is produced by interactions between 
these local processes and the large scale motion Held, and is 
not essential for the existence of the observed oscillation. We 
show that the interactions between radiative and convective 
processes over the warm tropical oceans can lead to low 
frequency oscillations. 

2. THE MODELS 

Two different models have been used to describe, with 
different degrees of physical detail, the local interactions 
among radiation, boundary-layer processes, and convection. 

In a highly simplified model (thereafter Model 1) the 
radiation process is approximated by cooling to space 
(Rodgers and Walshaw, 1966) and convection is included 
through the adjustment scheme of Manabe et al. (1965). 

In a more sophisticated model (Model 2), which is 
essentially a single-column version of the Colorado State 
University (CSU) general circulation model (GCM) described 
by Randall et al. (1991) and references therein, we use a full 
parameterization of solar and infrared radiative transfer 
(Harshvardhan et al., 1987), and the Arakawa- Schubert 
cumulus parameterization (Arakawa and Schubert, 1974). 

A fixed sea surface temperature (SST) is imposed as 
the lower boundary condition in both models. The latent and 
sensible heat fluxes from the ocean surface are determined 


using the bulk aerodynamic formula with a prescribed surface 
wind speed. 

Both models reach aquasi-equiiibrium state after a few 
hundred simulated days. The results for the last 500 davs of 
2000-day integrations are presented in the following section. 

3. RESULTS 


Fig. 1 shows the variation of precipitation amount from 
Model 1, with SST*301.5 K and a surface wind speed of 5 m 

s* 1 . 
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Figure I: lime hisory of the precipitation rate from Model 


An oscillation with a period near 34 days is evident (Fig. 2). 
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Figure 2: The spectrum of the precipitation data shown in 

Fig. 1. 



The fluctuations of atmospheric temperature and moisture 
during the first 100 days of this 500 day segment are plotted in 
the panels of Fig. 3. 
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Figure 4: The rime-history of the precipitation rate (minus 

its rime average) from Model 2. 
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Figure 3: a) The fluctuations of atmospheric temperature 

during a 100 day segment of the simulation with 
Model 1. The contour interval is 0.1 K. b) The 
fluctuations of atmospheric moisture during a 
100 day segment of the simulation with Model. 

The contour interval is 0.1 g kg' 1 . 


The temperature variations in the middle and upper 
troposphere are out of phase with those in the lower levels; 
warming in the lower levels is accompanied by a cooling in the 
upper levels and vice versa. Intense convecnon occurs while 
the upper levels are warming/moistening and the lower levels 
pe cooung/drying. Weaker convection occurs while the upper 
levels are cooling/drying and the lower levels are warming 
moistening. The variations shown in Fig, 3 illustrate a thermal 
destabilization as moisture and heat accumulate in the lower 
levels and radiative cooling occurs aloft, and a thermal 
stabilization as convection transpons both moisture and 
sensible heat from lower levels to the upper levels. 

The results from Model 2 are shown in Figs. 4 - 6. 


Figure 5: The spectrum of the precipitation data shown in 

Fig. 4. 
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Figure 6: a) The fluctuations of atmospheric temperature 

during a 100 day segment of the simulation with 
Model 2. The contour interval is 0.1 K. b) The 
fluctuations of atmospheric moisture dunne a 
100 day segment of the simulation with Model 2. 

The contour interval is 0. 1 g kg* 1 . 
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[n these runs we used an SST of 300.0 K and a surface wind 
speed of 5 m s' 1 . Model 2 produces oscillations of temperature 
and moisture that are quite similar to those of Model l. This is 
remarkable given the drastic differences in the formulations of 
the two models. As seen in the figures, sensible heat and 
moisture accumulate in the lower levels; this is followed by an 
episode in which moisture is transported upward by more 
intense convection. Again a cycle of thermal destabilization 
and stabilization occurs with a period of 53 days. 

Deep convection intensifies when moistening is a 
maximum, the upper levels are cool, and the lower levels are 
warm. Through both radiative and latent heating, more intense 
convection leads to a warming of the upper troposphere. It also 
moistens the upper troposphere through cumulus detrainment. 
The warming of the upper levels represents a thermal 
stabilization of the column, which suppresses convection and 
allows an accumulation of sensible heat and moisture in the 
lower levels. This accumulation is followed by a new round of 
more intense convection. 

Similar low frequency oscillations occur, with varying 
intensities, for a wide range of SSTs and surface wind speeds. 
More intense oscillations are associated with wanner SSTs. 

4. DISCUSSION AND CONCLUSIONS 

Our strategy now is to simplify Model 2, in a step-by- 
step process, to identify which physical processes are 
important for the oscillation and which are not. The fact that 
the oscillation does occur in the much simpler Model 1 can 
guide this process. We have already determined, through 
experiments with Model 2, that the process does not depend in 
any important way on the radiative effects of the clouas or on 
the variability of the planetary boundary layer depth. 
Additional experiments are in progress. 

We offer the following tentative explanation for our 
results: Conditional instability can be generated, in our l-D 
models, by radiation and boundary-layer processes. Both 
moist convective adjustment and the Antkawa-Schubert 
scheme instantaneously release conditional instability, i.e., 
they force the models to occupy moist neutral states. In both 
models, the existence of moist instability depends on the 
temperature lapse rate and also on the amount and distribution 
of water vapor. The oscillations produced by the models may 
be due to the existence of two '‘attracting” moist-neutral states: 
the first with a stronger temperature lapse rate and a lower 
moisture content, and the second with a weaker lapse rate and 
more moisture. The convective and radiative tendencies 
associated with the first state drive the system towards the 
second, and vice versa. The period of the oscillation is that 
associated with the radiative cooling, and perhaps also by the 
time needed for moisture to accumulate in the column. 

Our results show that a low frequency oscillation can 
be excited through radiative-convective interactions in the 
presence of energy exchange at the sea surface. This 
oscillation may well be responsible for the observed stationary 
low frequency fluctuations in the tropics. It may also provide 
a foundation for studies in which a pre-specified oscillatory 
forcing is driving the tropical dynamic process (e.g.. Gill, 
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1. INTRODUCTION 

Since better understanding of 30-60 day oscillation can 
lead to improvements in medium and long range weather pre- 
diction, this phenomenon has attracted considerable attention. 
The 30-60 day oscillation was identified by Madden and Julian 
(1971, 1972) as the eastward propagation of a wavenumber 1 
disturbance in tropical, tropospheric zonal wind and surface 
pressure fields. The range of the period of this disturbance 
corresponds to a zonal phase speed of 10 to 15 ms' 1 . Other 
observational studies have documented disturbances with a 
similar period of 30 to 60 days in other physical fields. These 
findings include the oscillations in atmospheric angular mo- 
mentum (Anderson and Rosen, 1983, and Risbey and Stone, 
1988), convective activities and cloudiness (Yasunari, 1980), 
outgoing long wave radiation over the Indian- Western Pacific 
Oceans (Murakami et of,, 1986, among others) and perturbs 
tions in precipitation amount over the same area (Hartmann 
and Gross, 1988). In addition to atmospheric phenomena, 
30-60 day fluctuations of sea level height (Enfield, 1987, and 
Mitchum and Lukas, 1987) and sea surface temperature (SST) 
(Enfield, 1987) have also been found. From these studies it is 
obvious that the 30-60 day oscillation is a complex process 
involving disturbances in both the tropical atmosphere and 
ocean. 

Despite the multitude of observational studies, a satis- 
factory theoretical explanation of the 30-60 day oscillation has 
yet to be agreed upon. Since Gill's (1982) finding, that the 
phase speed of the tropical waves can be significantly reduced 
due to reduction of the static stability (hence the equivalent 
depth) of the atmosphere by latent heat release, the linear 
tropical wave theory (Mateuno, 1966) after being refined by 
the diabatic heating parameterization, has been widely ap- 
plied in explanation of the 30-60 day oecillation. For example, 
Lau and Peng (1987) proposed the ‘mobile wave CISK’ mech-’ 
aniam and Chang and Lim (1988) developed the ‘Kelvin wave 
CISK* theory. Although the selected or the newly generated 
waves in these theories were indeed amplifying as propagating 
eastward the phase speed of these waves was crucially depen- 
dent on the diabatic heating profile used in these models. In 
order to have a phase speed close to the observed one the heat- 
ing profile has to have a peak in the lower troposphere. This 
heating profile, however, is more representative of situations 
for the mid-latitude than of the tropics (Reihl and Malkus, 
1958). 

Attempts were also made to find new modes which 
could be responsible for the observed 30-60 day oscillation. 


Anderson and Stevens (1987), for example, specified several 
slow poleward propagating unstable modes in their 2-D eigen- 
value model which included an equatorial symmetric Hadley 
circulation. 

Studies in which the intrasessonal oscillation was re- 
garded as an atmospheric response to a sporadic external heat- 
ing with a 30-60 day period is another category in the ap- 
proach. By postulating an oscillatory heating as the forcing, 
many experiments with different atmospheric models were able 
to produce the expected low frequency oscillatory features ( Yar 
magada and Hayashi, 1984, and Salby and Gacia, 1987). The 
lack of an explicit instability mechanism in this category of 
explanation remains questionable, however, if one views the 
30-60 day oscillation as purely an intrinsic atmospheric oscil- 
lation. 

The sporadic ‘external’ forcing may turn into an ‘in- 
ternal’ one when the ocean and atmosphere are considered 
as one system. Krishnamurti et si. (1988) has found the 
observational evidence of the air-sea interaction on the time 
scale of 30-60 days. Emanuel (1987) and Neeling et a t (1987) 
looked at the evaporation-wind feedback effect in an ocean- 
atmosphere system. They showed that the low frequency os- 
cillation can be amplified by drawing energy from the warm 
ocean surface through this feedback. Although the existence 
of the oecillation does not depend on the evaporation-wind 
feedback (Neeling it «/., 1987) their studies certainly provided 
a different perspective for understanding this problem. 

In this study, we propose a new theory explaining the 
origin of low frequency oscillations in the tropics. In our mod- 
els both the lower atmosphere and upper ocean layers are in- 
cluded and are treated as a fully coupled system. Based on 
the model results, we suggest that the observed 30-60 day os- 
cillations in the tropical tropoephere and in the ocean are the 
manifestation of the cloud-radiation, atmosphere-ocean feed- 
back processes, and that their realization depends on the hy- 
drologic cycle in this system. 

In the next section, the basic physical mechanism for 
the 30-60 day oscillation in terms of the hydrologic cycle is 
presented. Sections 3, 4 and 5 describe three different models 
and their results. Our conclusions are summarized in Section 
6 . 

2. MECHANISM 

The tropical atmosphere above the warm ocean surface 
is heated through the energy exchanges with the ocean. Mean- 
while evaporation from the ocean moistens the lower part of 



the tropical troposphere. The warmer the ocean mixed layer, 
the more heat and more moisture are added into this part of 
the atmosphere. 

The removal of the accumulated ‘extra* heat and mois- 
ture from the tropical atmosphere occurs through deep cu- 
mulus convection (Riehl and Malkus, 1958). This convection 
is crucial in the hydrologic cycle of the tropical atmosphere- 
ocean system. When explaining the conditional instability in 
the tropical atmosphere, Ooyama (1969) showed that the verti- 
cal temperature and moisture structures in the tropics possess 
a conditionally unstable character before the onset of convec- 
tion. However, this instability is realized only when the lower 
troposphere is sufficiently moistened, through surface evapora- 
tion and/or large-scale convergence. This characteristic of the 
tropical atmosphere is important as it highlights the key role 
played by the moisture in the occurrence of tropical convection. 
In our conceptual models, presented in the next two sections, 
we will set a critical moisture value as the model atmospheric 
moisture tolerance’ to represent this thermodynamic feature 
of the tropical atmosphere. This treatment principally follows 
the philosophy in parameterizing cumulus convection (Kuo, 
1975, and Stevens and Lindzen, 1978) and is based on empiri- 
cal evidence such as the observed moisture changes in synoptic 
scales wave disturbances (Reed and Recker, 1971). 

When the conditional instability is being released, the 
convection not only diminishes the extra moisture in the at- 
mosphere through precipitation but, through cloud-radiation 
interaction, also reduces the sea surface temperature by atten- 
uating the solar radiation energy available at the ocean surface. 
The result of this is to return the extra water substance in the 
atmosphere back to the ocean and to reduce the sea surface 
temperature to its ‘climatological* value. 

In short, a natural cycle represented by the hydrologic 
process in the atmosphere-ocean system is spontaneously car- 
ried on in the tropics. Our model shows that, simply restricted 
by the ‘moisture tolerance’, such a cycle possesses an intrasea r 
sonal time scale. 


mixed layer through evaporation and the removal of this con- 
tent by precipitation; horizontal convergence and divergence 
are neglected. 

The temperature of the atmospheric mixed layer is 
characterized by a near-surface atmospheric temperature T a = 
T t - AT(°C). Precipitation is assumed to take place when 
the ‘atmospheric moisture tolerance’ is achieved. As an initial 
approach we specify this critical condition as that when the 
mixing ratio of the atmospheric mixed layer reaches 90% of 
the saturation value at its temperature, TV 

The amount of moisture evaporated in each time in- 
terval, At, is calculated by 

AE = f E{t)dt, (1) 

Jtkt 

where, 

£(0 = P.CoMlq^T.) - q atm (t)\ (2) 

(Neeling tt a l 1987). In (2), p* is the mean density of the 
atmospheric mixed layer, C D is the surface drag coefficient and 
is treated as a constant of 0.001. The parameter u, is the low 
level wind at 10 meters above the sea surface, 9 Aa *(7*) is the 
saturation mixing ratio of the atmosphere at temperature T, 
and <h tm (t) is the actual mixing ratio of atmospheric mixed 
layer at time t. The value of <j a tm is updated at time step 
t + At by an increment A q given by 

<•> 

Here P represents the precipitation rate and is assumed as a 
constant, 2.5mm/day (Wang, 1988) during the moder§ rainy 
period and H is the mixed layer depth. 

By using a mean atmospheric moist state and the pre- 
cipitation rate, the model is integrated with a one hour time 
interval for 1000 days. 

b) Model Result and Sensitivity 


3. A SIMPLE MODEL WITH FIXED SST 
a) The Model 

In this section we demonstrate the existence of the hy- 
drologic cycle in a simple atmosphere-ocean system with fixed 
ocean parameters. The schematic structure of this model is 
given in Fig. 1. As shown in Fig. 1 only thermodynamic pro- 
cess are considered in this model. Our attention is to describe 
the accumulation of the extra water content in the atmospheric 
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Fig. 2 shows the results of two model integrations. 
Model parameters used in these integrations are listed in the 
figure caption. The displayed field, extra water content, in- 
dicates the column moisture amount above an assumed back- 
ground state, which corresponds to a mixing ratio (? a tm) of 
15 gjkg in the atmospheric mixed layer. 

An ad hoc parameterization of the low-level wind, in- 
tended to represent the higher winds and associated transfers 
between atmosphere and ocean during convection, 



Fig. 2: Predicted hydrologic circulation in the model atmosphere- 
ocean system. The model low level wind, u # , used in (a) is 
Fig. 1: Schematic structure of the model (MD1) with a a constant of 2.0 m/s, and in (b) is that given in Eq. f4), 

fixed sea surface temperature (SST). with u* = 2 m/s, u m = 12 m/s, and At* = 13 days. Other 
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is used in getting the result in Fig. 2b. In this expression u fc 
is the mean background wind speed and t e is the time when 
the atmospheric moisture tolerance is achieved. The maximum 
wind speed, u*, is assumed to occur during the prescribed time 
interval At m when deep cumulus convection is present. This 
time increment is prescribed to be 2 days, 

Ftom the results shown in Fig. 2 the existence of the 
hydrologic cycle operating in this simple atmosphere-ocean 
system is clearly evident. Through the processes of evapora- 
tion and precipitation, which are constrained by the moisture 
tolerance and its implied conditional instability, the hydro- 
logic cycle is naturally introduced in this system. Further- 
more, through a comparison of these two results, we see that 
the presence of this cycle is independent of the detailed struc- 
ture of the low level wind. This independence suggests that 
the oscillation in the low level wind field may be considered a 
consequence of the hydrologic cycle. Since the time evolution 
of the heating profile in the tropical troposphere is directly re- 
lated to the different phases of the hydrologic cycle, this cycle 
provides an external heat/energy source for the atmoephere. 
This source disturbs the energy state of the tropical atmo- 
sphere to which the wind and mass fields of the atmosphere 
then adjust. Thus this hydrologic cycle is a plausible mech- 
anism for initiating oscillations of similar frequency in wind 
and mass fields in the tropical atmosphere. The influence of 
the varying low level wind field on the frequency of the hydro- 
logic cycle is also noticeable in Fig. 2, although it is does not 
determine the existence of the cycle per it. 

Results of model integrations with different model pa- 
rameters are presented in Fig. 3. These results were obtained 
with the same low level winds as used to produce Fig. 2b. A 
significant characteristic shown in this figure, is that the period 
of the resultant hydrologic cycle in this system increases with 
increasing model SST. This result is consistent with the obser- 
vation that low frequency oscillations are preferred in regions 
of warmest SSTs (Webster, 1987). 

The model sensitivities to different low level wind 
speeds and values of precipitation rate are also illustrated in 
Fig 3. 

4. COUPLED INTERACTIVE ATMOSPHERE-OCEAN 
MIXED LAYER MODEL 

The existence of the hydrologic cycle and its relation- 
ship to the low frequency oscillation in the tropical troposphere 
have been shown by the one-way interaction model (hereafter 
MD1) of Section 3. Since a fixed SST was applied in MD1, the 
question of how this cycle as well as its time scale is influenced 
by a time varying SST still remains. In this section we address 
some aspects of this issue. 

a) The Model 

The model structure and the m^jor physical processes 
are schematically illustrated in Fig. 4. As in MD1, we focus 
on the thermodynamics of this atmosphere-ocean system and 
ignore any detailed dynamics. The model's ocean mixed layer 
process follows that in Kraus and Turner (1967), and Denman 
(1973) 

In this model to incorporate in the model SSTs which 
vary in time, a detailed energy budget in the ocean mixed 
layer is required. As the only energy source in this budget, 




Fig. 3: Predicted period of hydrologic circulation with dif- 
ferent model conditions. The abscissa shows the 
temperature difference between the ocean and the 
atmosphere, (a) and (b) show the model sensitiv- 
ities to the changes in precipitation rate P [unit: 
mm/dayj and the model low level wind u, unit: 
m/sl, respectively. In (c) and (d), the 3 groups of 
results, e.g. A*, 2?< and C, are for 3 different P’s 
used in the model, i.e. 2.5, 3.0 and 3.5 mm/day. 
The subscript », from 1 to 3, stands for different 
low level wind conditions of 2.0, 3.0 and 4 0 m/s 
Lines are not plotted in (c) and (d) when a bal- 
ance between E and P is the case and therefore 
no osculation occurs. 

the input solar energy at the ocean surface needs to be in- 
cluded. This input further varies with presence or absence 
of cumulus convections through cloud-solar radiation interac- 
tion. On the other hand, the energy losses through infrared 
radiation (IR) from ocean surface, and sensible and latent heat 
exchanges with ths atmosphere are also considered. Since ex- 
plicit prediction of key processes, such as cloud generation and 
interaction of clouds with solar radiation, is beyond the scope 
of this mechanistic study, we parameterize their effects. To ex- 
pttm the cloud influence on the solar radiation, for example, 
we assume different values of the solar irradiance at the sur- 
face under varying convective cloud cover conditions. These 
values are obtained from observations and theoretical studies. 
A similar method is also applied to represent the IR cooling 
process. The main reason for neglecting the detailed treatment 
of these processes is to avoid ths uncertainties associated with 
complex processes involving multiple scales, toms of which are 
poorly understood, and therefore to concentrate more on the 
basic physical mechanisms. 







Fig. 4: Schematic structure of the fully coupled atmosphere- 
ocean mixed layer model (MD2). 


In the model atmosphere, the process is similar to 
MD1. The moisture amount transferred into the model at- 
mospheric mixed layer is calculated again using the bulk aero- 
dynamic formula 

£(0 = P«c D MOI (*.i [T.[t)\ - *•«(*)) • (5) 

Here, q,at[T,(t)\ is time dependent as SST varies with time in 
this model. A Bowen ratio of 0.2 is used to obtain the value of 
the sensible heat flux required as an input to the ocean model. 

The solar irradiance at the ocean surface are specified 
by the step function 

p _ f R ?w » non-convective; 

~ \ R c 0Wt convective. 1 ' 

The values in (6) are based upon the satellite and observational 
data analyses by Gautier and Katsaros (1984), Gautier (1986) 
and Gautier (1988). With these values, the bulk features of 
the influence of deep cumulus convection on the solar radiation 
transfer in the atmosphere have been included. Likewise the 
step function for IR change is 


p _ / K?/i> non-convective; 

1 R \ RJ/i, convective, ' ' 

where the values for R* n and R r in are based on a study by 
Stephens et al. (1981). 

b) Model Results 

Table 1 lists some of the model experiments we con- 
ducted. The physical parameters, their values, and the model 
predicted quantities are included in this table. 


Since only thermodynamic and turbulent mixing pro- 
cesses in the two mixed layers are considered, only a few 
of these experiments reach a steady or an equilibrium state 
through the period of integration. Fig. 5 shows the result from 
one of the experiments (experiment 1 in Table 1) in which an 
equilibrium state has been reached. From this figure we see 
that the oscillation of the extra water content in the atmo- 
spheric mixed layer still experienced when the model has vary- 
ing SST. This oscillation, through the coupling and feedback 
with the solar radiation and hence the surface energy process, 
influences the processes in the ocean mixed layer (Figs. 5b, c 
and d). The time scale of these oscillations predicted by the 
model is 33 days. 

As we have discussed in Section 2, in the process of 
the hydrologic circulation a pulsating energy disturbance is re- 
leased in the convection phase. This disturbance then causes 
dynamic perturbations in the atmospheric momentum and 
mass fields (Salby and Garcia, 1987). The observed oscillation 
in the atmospheric relative angular momentum (Anderson and 
Rosen, 1983) may also be related with this pulsating energy 
disturbance. On the other hand, this energy disturbance, in 
addition to the change in the solar energy input, affects the 
energy budget of the ocean mixed layer. As a consequence 
of the hydrologic circulation in the atmosphere-ocean system, 
the 30*60 day oscillations therefore appear in both the tropical 
atmosphere and ocean. 

Because the thermodynamic processes depend on SST 
in a highly nonlinear fashion, the results illustrated in Table 1 
differ from those obtained from MD1. In Fig. 3, a higher SST 
in the one-way model causes a dramatic increase in the period 
of the model’s hydrologic circulation. Results from the fully 
coupled model, however, are not as simple. Using 29 °C as the 
initial SST, the hydrologic cycle reaches a steady period of 35 
days (experiment 2 in Table 1). This period is quite close to 
that predicted in experiment 1 where the initial SST is 28°C. 
However when a colder SST is used to start the model run, 
the predicted result are distinctly different (Fig. 6). 

Results presented in Table l also show the model sen- 
sitivity to other model parameters. 

5. ONE-DIMENSIONAL GCM AND RESULTS 

To test our hypothesis in a more complete manner dif- 
ferent experiments have also been conducted with the one- 
dimensional version of the UCLA/CSU GCM (Randall and 
Dazlich, 1989). This 1-D model includes the full radiation and 
moist physics parameterization of the GCM, and is coupled 
to a slab of ocean of fixed depth. The slab ocean does not 
exchange energy with deeper ocean layers; its temperature is 
controlled entirely by the surface energy flux. 

At this stage some very preliminary results from ex- 
periments with this 1-D GCM model have been obtained and 
are shown in Figs. 7 and 8 . These results support the ma- 
jor findings from our previous conceptual models. When run 
with a ocean slab depth of 60m, the model produces obvious 
oscillation of sea surface temperature, with amplitude 0.4°C 
and a period of 60 days (Figs. 7a and 8a). This oscillation is 
accompanied by 60-day fluctuations of the surface energy flux, 
precipitation rate and cloudiness (Figs. 7b, 7c and 8b). The 
changes in the surface energy flux are due to changes in the ab- 
sorbed solar radiation, which are controlled by the cloudiness 
fluctuations or model convections. 

During the warming phase of the oscillation, the sea 
surface temperature and the convective precipitation rate 
gradually increase. Eventually the cloudiness increases to the 
point that the solar radiation absorbed by the ocean is reduced, 
allowing the ocean to begin cooling. This cooling is followed 
by a reduction in the model convective activity, which leads 
to a reduction in the cloudiness. The absorbed solar radiation 
then increases, and the model hydrologic cycle is re-initiated. 



TabU 1< Modal Experiments and Results - Among the output quantities, r is the 
predicted period of the model hydrological circulation. Zq is the maximum extra water 
amount accumulated in the model atmosphere before convection starts. AT and AA are 
respectively the amplitudes of changes in the SST and depth of ocean mixed layer in a 
hydrological circulation. When a steady state is not reached in a model run, these quantities 
are given from a nominal cycle in model integration. 


Physical 

I Experimental Value | 

Parameter 

1 

2 

3 

4 

5 

6 


250 

250 

250 

250 

250 

260* 


120 

120 

120 

120 

120 

105* 


40 

40 

40 

40 

40 

40 

RiSfiSi 

40 

40 

40 

40 

40 

40 

Klilifffi! 

I9BSI 

mm 

m 

m 


(2, 12, 13.0) 

P (mm/ day) 

2.5 

3.0* 

2.5 

2.5 

2.5 

2.5 

EHS83 

28.0 

29.0* 

26.0* 

28.0 

28.0 

28.0 

Emu 

25.0 

26.0* 

23.0* 

25.0 

25.0 

25.0 




Model Predicted Result 



r (d*y.) 

55 

35 

tha 5th: 25 
tha 9th: IS 

tha 5th: 33 
tha 20th: 37 

tha 5th: 33 
tha 20th: 26 

31 


1.49 

1.51 

tha 9th: 0.35 

tha 20th: 1.70 

tha 20th: 0.80 

1.40 


1.68 

1.63 

tha 9th: 0.53 

tha 20th: 1.67 

tha 20th: 1.41 

1.83 

1 

21.85 

21.39 

tha 9th: 13.51 

tha 20th: 21.48 

tha 20th: 19.50 

24.93 


*: Parameter value different from that in the controlled experiment (experiment 1). 
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Fig. 5: Model results from experiment 1 in Table 1. 
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Fig. 6: Temporal variation of the extra water content in 
the model atmosphere from experiment 3 in Table 
1 . 
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Fig. 7: Temporal variations of (a) SST, (b) solar radia- 
tion absorbed at the ocean surface, and (c) daily 
precipitation from the 1-D GCM. 



6. CONCLUDING REMARKS 
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Fig. 8: Square of the amplitudes of the Fourier harmonics 
of (a) model SST, and (b) model precipitation. 
Peaks at about 60 days are outstanding. 


In this study we have shown that the 30-60 day oscilla- 
tions in the tropical atmosphere and ocean can result from the 
hydrologic circulation in a coupled atmosphere-ocean system. 
This circulation is intrinsic to this system and is a consequence 
of the energy exchange and energy balance in this system. 

The time scale of this circulation can be determined 
in a simple coupled atmosphere-ocean model when the role 
of moisture played in the atmosphere is incorporated in the 
model. By parameterizing this role in terms of the ‘model 
atmospheric moisture tolerance 1 in our conceptual models we 
showed that the hydrologic cycle is spontaneously carried on 
in the tropica) atmosphere-ocean system. Its existence was 
basically independent of the variations of the low level wind. 

Our hypothesized relationship between the hydrologic 
circulation and the intraseasonal oscillation in the tropical at- 
mosphere and ocean was also examined in a 1-D GCM. The 
model results were positive in supporting this relationship. 
Meanwhile, the analysis of the thermodynamic process in the 
1-D GCM illustrated the mechanism we proposed. 

Based on our results we conclude that the hydrologic 
cycle in the tropica) atmosphere-ocean system is fundamen- 
tal to the observed intraseasonal oscillations in the tropics. 
The physical restriction of the accumulation of latent heat 
in the lower tropica) troposphere and the reduction of the 
source supply of energy to the tropical atmosphere-ocean sys- 
tem, through cloud-solar radiation interaction, are crucially 
important to the presence and time scale of this hydrologic 
cycle. 
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SUMMARY 


Research Triangle Institute (RTI) is a subcontractor to Colorado State University 
(CSU) on the National Aeronautics and Space Administration (NASA) project "Linear and 
Nonlinear Aspects of the Tropical 30-60 Day Oscillation: A Modelling Study," contract 
number NAG5-1 122. The proposed three year program started in March 1989 and will 
conclude in February 1992. In this integrated program, RTI has to develop a 5-level 
general circulation model (year 1 and 2) and design numerical experiments to analyze the 
model simulation of the 30-60 day oscillation (year 3). During the second year of the 
project (March 1990 - February 1991) RTI completed the development of the 5-level 
general circulation model. Specific tasks that have been accomplished during this year are: 

1 . Test the time difference method used for model integration; 

2 . Design the radiation scheme that is used in the model; 

3 . Develop test runs for the model simulation for the month of January (Northern 
Hemisphere winter). 

RTI principal investigator, Dr. Adel Hanna has maintained full interaction with the 
co-principal investigators. Dr. Graeme Stephens, CSU and Dr. Duane Stevens, University 
of Hawaii (UH). Also, he has established scientific cooperation with Dr. James Hack of 
the National Center for Atmospheric Research (NCAR), who is responsible for developing 
the next generation of the Community Climate Model (CCM). 
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1.0 INTRODUCTION 

Modelling analyses provide an appropriate tool to understand atmospheric circulation 
especially on a global scale where the observational network is not dense enough to account for all 
the scales of atmospheric interactions. The intraseasonal (30-60 day) oscillation has received a 
great deal of attention due to its possible impact on global climate change as an intrinsic mode of 
the atmospheric interannual variability. Both observational and modelling studies of the 
intraseasonal oscillation have proven to be useful in identifying the origin of its structure. The 
oscillation is mainly driven by nonlinear wave interaction of the Rossby and Kelvin waves in the 
atmosphere. These waves are excited by the latent heat release associated with intense convective 
systems over the tropical oceans. On the other hand the latter is influenced by the sea-surface 
temperature pattern. 

This research project is an integrated program to study the linear and nonlinear aspects of 
the 30-60 day oscillation. The nonlinear dynamics are investigated through the analysis of a five 
level primitive equation global climate model. The model is of medium complexity structure to 
understand and isolate specific components of the atmospheric circulation. However it includes the 
forcing mechanisms and physics to stimulate the large scale atmospheric circulation and seasonal 
cycle through interactions of orographic forcing, radiative heating, convective and nonconvective 
release of latent heat. The development of the model in both the finite difference and spectral forms 
gives the opportunity to study the interactions between different wave numbers and also the effect 
of spectral truncation on model simulation. For example, one of the questions of basic interest to 
this study is the relationship between zonal asymmetries (especially wave number 1) and zonally 
symmetric mode (wave number 0) in producing energetics necessary for maintain ing the 
oscillation. 

This report represents a description of the 5-level general circulation model. It also 
includes a description of the projected experiment using the model. 

2.0 SYSTEM OF EQUATIONS 

The 5-level primitive equation model has been developed. The atmospheric model is global 
using sigma (= p/p,) as the vertical coordinate. The horizontal model resolution is spectral using 
spherical harmonics rhomboidal truncation at wave number 15. The prognostic variables are the 
stream function (¥), the velocity potential (X), surface pressure (p,), temperature (T), and the 
mixing ratio (q). The stream function (\j f) and velocity potential (x) are expressed in Laplacian 
forms (vorticity- and divergence - 5). Diagnostic variables are the geopotential height (<I>), and 

the vertical velocity (o). The following system of equations represents the basic equations for the 
5-level atmospheric model (see for example, Bourke, 1974; Haltiner and Williams, 1980). 
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The zonal velocity, u, and the meridional velocity, v, are calculated using the stream 
function and velocity potential. The diabatic heating rate, Q , and moisture source, M, are 
calculated using surface (land and ocean) and atmospheric fluxes. Other parameters in the above 
system of equations are: E - the Kinetic energy; f - the coriolis parameter; <t> - the latitude angle; Cp 
- the specific heat; R - the gas constant; x x and ty - the horizontal stresses; g - the acceleration of 
gravity; a - the radius of the earth; and Di - the diffusion of appropriate parameters. The overbar 
presents the vertical integration through the depth of the model atmosphere, or from the top of the 
model atmosphere (<r=0) to surface 0=1. The prime presents the deviation from the horizontal 
average defined by the subscript 0. The temperature (T) and geopotential (<J>') fields are expressed 
in the form of a departure from the horizontal average in order to apply a semi-implicit time- 
differencing scheme. 

3.0 HORIZONTAL TRUNCATION; SPECTRA L METHOD 

The horizontal representation of the dependent variables is in terms of a truncated series 
using spherical harmonics as basis functions. Nonlinear terms are transformed to the spectral 
space using the full transform method. The variables are expanded in the form 
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where ( ) n m identify the harmonic coefficients, ^ is longitude, M- is the sine of latitude, m is the 
zonal wave number, n - 1 m 1 is a meridional wave number, M is the highest zonal wave number 
retained in the truncated series, and J is the highest value of n - 1 m 1 retained in the truncated 
series. Y B m are spherical harmonic function defined by 


Y^ = P? 00 e*"*. 


3.12 


where P™ are the associated legendre functions of the first kind. 

The horizontal truncation using the spectral method gives the model flexible efficiency in 
terms of using different zonal wave numbers to test different horizontal resolutions and the 
interaction between the basic flow and certain wave numbers. In the current program we assume 
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rhomboidal truncation at wave number 15, with 40 Gaussian latitudes and 64 longitude grid points 
at each latitude. 

The rotated zonal wind U and meridional wind V components, are diagnosed from the 
stream function and velocity potential harmonics using the forms; 
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4.0 VERTICAL FINITE DIFFERENCE AND I NTEGRATION 

In the vertical, the variables are staggered in a in such a way to calculate vertical derivatives 

and the vertical velocity, & between layers (Figure 1). Levels 1 to 5 are located at layer midpoints. 


Vertical advection terms of momentum, temperature, and moisture at level k is calculated using the 
form 
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The vertical velocity is zero at the top and surface levels of the model's atmosphere. 
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In Finite difference form, the vertical velocity a (equation 2.6) is given by 
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where K is the number of levels (K=5), and 

Sj=(* 2 x). 45 

5 0 EA RAMETFiRIZATION OF PHYSTrAT. p r ocesses 

The diabatic heating rate.Q, is a composite of the radiative heating, latent heating and 

surface heat fluxes as boundary interchange between the surface (ocean and land) and the 
atmosphere. 

The surface heat and moisture fluxes are calculated using a bulk aerodynamic formula. The 
surface sensible heat flux (QJ is given by: 

Q* = P*CpC d |V s |(T g -T 5 ) 51 

where T 5 is the air temperature at level 5, p, is the density of air, I V 5 I is the windspeed, and C d is 
the drag coefficient Ground temperature T g over land is determined from the surface energy 
balance, whereas on the ocean T g is specified from climatological values. 

The evaporation from the surface to the lower layer of the atmospheric model is given by: 

^“PiCdl V 3 IGW[q 4 (T g )-h s q J (T 5 )], 5 . 2 

where h s is a measure of the relative humidity in the lower layer, q s is the saturation mixing ratio 
and GW is a wetness parameter that is equal to unity on water surface. 

The amount of latent heat lost from the surface due to evaporation is given by: 

Qe= , 5 3 

where L is the latent heat of condensation. 

5,1 L ARGE-SCALE CONDENSATIO N AND CONVECTIVE 
adjustment 

The convective adjustment method used in the model follows the algorithm 
developed by Manabe et al., (1965). Each level is tested for supersaturation and hydrostatic 
stability. Excess moisture in layers which are stable and supersaturated is condensed as large-scale 
precipitation with release of latent heat in the layer. Supersaturated and convectively unstable 


6 



layers are adjusted in such a way to conserve moist static energy through the release of latent heat 
and the vertical transfer of heat and moisture to the successive layers aloft The procedure and tests 
for convective adjustments are explained in the following; 

a) Large-scale condensation only 

r£ 1 < p v 

7) n ^ r s *m 

/ \ < * 
q+5q = q, (T + 5T, a) 

Cp 8T + L8q = 0 

b) Moist convection and large-scale condensation 

r£ 1 

!|.(T + 5T,q + 5q,<y) = 0 , 
q + 5q = q, (T + 5T, or) , 

J (cp5T + L8q) da = 0 
o 

where r is the relative humidity, Y m is the moist adiabatic lapse rate, 8q, and 8T are the adjustments 
to water vapor mixing ratio, and temperature respectively. 

The atmospheric moisture, in the model, is described in terms of the water vapor mixing 
ratio, q, (Equation 2.5). In physical space, this quantity is non-negative by definition. However, 
when using the spectral method the finite truncation in wave-number space generates negative 
values when transformed back to physical space. This misrepresentation is known in Fourier 
analysis as Gibbs phenomenon. Various formulations of the moisture equation using other 
required positive variables, such as the relative humidity, are also having the same problem of 
transformed negative values in the physical space. To overcome this problem, some spectral 
models are using the vertical borrowing method to augment the negative values to zeros. 
However, if there is no available moisture in the upper levels of the model, these negative values 
are simply increased to zeros. This procedure causes a spurious source of water vapor that affects 
the model's global moisture budget Royer (1986) has suggested a method based on horizontal, 
borrowing in such a way to maintain the conservation of the global moisture. The horizontal and 
vertical moisture borrowing methods are used in the 5-level model simulation. 
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5.2 SOLAR RADIATION 

The calculation of radiative energy of the model is similar to the one described by 
Hanna and Stevens (1984). The incoming solar radiation at the top of the model's atmosphere is 
calculated as a function of daily mean zenith angle (Wetherald and Manabe, 1972). The albedo of 
the atmosphere r t , is calculated, taking into consideration the observed zonal mean amount of 

clouds. 

The solar radiation absorbed by the atmosphere, S r , is given by: 

s,.xd-Os, 58 

where x is the absorptivity of the atmosphere, and S is the incoming of solar radiation at the top of 
the atmosphere. 

The net solar radiation absorbed by the earth's surface (land or ocean) is given by: 

S, = (1 • x ) (l * r a) (l • rj s, 5 9 

where r, is the albedo of the earth’s surface. Over the oceans, r, is assigned to .1, while over the 
continent, it is assumed to be a value between .1 and .8, depending on the type of surface (Gates, 

1975). 

5.3 LONG-WAVE RADIATION 

The calculation scheme of the longwave radiation cooling rates of the atmosphere 
uses empirical longwave radiative emissivities calculated from climatological monthly mean 
humidities, cloudiness, and cooling rates and the model simulated temperature (Corby et al., 

1977). 

The net flux of longwave radiation in the atmosphere is given by: 

= -eBT 4 , 5.10 

' 'Lw 

where c is the heat capacity of the atmospheric layer, B is the Stefan Boltzmann constant and e is the 
longwave emissivity. 

The longwave emissivities at each level are calculated using climatological zonal mean 8 50- mb 
relative humidities for the simulated month (season). The zonal mean 850-mb relative humidities (Qsso) 
are specified at Gaussian latitudes. Empirical relations are used to calculate the longwave emissivities at 
different levels. The following represent the emissivities, at the model five levels, used for the January 
simulation: 


3T 
c dt 
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At the surface. 


Ej = .132 

e 2 = .17 + .0075 q 830 (<J>), 

£ 3 = .171 + .0027 q 830 (<», 
e 4 = .198 - .0042 q 83Q (#, 
e 5 = .147 -.0032 q 85Q (#. 

e g (<t>) = .4 - .08 (q^o) t 1 ’ c i) , 


where Cj is the cloud amount taken to be *» .5. The sum of effective emissivities for the model's 
atmosphere is given by 

S(« = -974 - .004 q g50 (#. 


A correction factor is applied to coefficients Ck (k = 1 to 5) so that 


t 

k» 1 


e k ((» = S«»-e g «». 


The net flux of longwave radiation at the earth’s surface is given by: 


L 0 = -e g BT g 4 

where e g is the surface longwave emissivity. Over open ocean, the surface temperature is 
specified as climatological values. Over land, snow, and sea ice, the surface energy is computed 
assuming thermal surface-energy balance (Holloway and Manabe, 1971). 


" ®g ® + S* - Qj - Q« + Qj^ - 0 

The term Q ice represents the heat conduction into and out of pack ice and is defined as 


5.12 





(VV) ovarpckre \ 
ever 41 exher surfaccy 


5.13 


where c l is the thermal conductivity of ice (* 2.092 wm’ 1 k* 1 ], I is the ice thickness (=*2m), and T gw 
= 271 .2°K is the temperature of sea water below the pack ice. The implicit equation (5. 12) for T g 
is solved by the Newton-Raphson interactive procedure (see Hanna and Stevens, 1984). 
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6.0 TTME DIFFERENCE METHOD 

The system of primitive equations that describes the model include the fast-moving gravity 
waves; which imposes a certain limitadon on the time step used in numerical integration to maintain 
computational stability if a centerd time method is used. For this reason the fast moving gravity 
waves are treated semi-imp licitly. Assume a general prognostic equation of the form 

^.,F„(Z)+r(Z) + F (Z) 61 

9t G 5 

where F G is the vertical diffusion fluxes, 

F s is linear horizontal diffusion calculated in spectral space, 

T is all the remaining terms. 

The three step time method is applied using the following equations (Williamson et al. 1987) with 
each identified by its "TIME" stage in accordance with the flow diagram ( Figure 2). 


TIME 1: 

Z n+ 1 = Z ' 1 + 2AtF 0 m 0 

6.2 

TIME 2: 

z ~n + i _ z ^. + 2At f(z ^ , Zf 1 , Z™, Z"’* 1 ) 

6.3 

TIME 3: 

Z*" 1 * 1 = Z**"* 1 + 2AtF s (Z**^ 1 ) 

6.4 

The bar defines a time filter 

Z n = z"+ a (z"' 1 - 2Z"+ z 1 * 1 ) 

6.5 


where oc = .06 

7.0 MODEL STRUCTURE AND PROGRAM DESCRIPTION: 

The Fortran code of the model consists of about 4,000 statements that execute forty-two 
functions and subroutines. The code is designed to run on the VAX/VMS at RTI. However, the 
code is designed also to be used on the CRAY computer series. The simulation is designed for the 
month of January (Northern Hemisphere winter) where the climatological values of sea-surface 
temperature and 8 50- mb relative humidities are used to calculate surface fluxes and long-wave 
radiation emissivities needed to integrate the model. We categorize the model subroutines into four 
groups: 
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1) . Preparation Subprograms; 

2) . Transform Subprograms; 

3) . Model Physics Subprograms; 

4) . Time Integration Su b progr a ms 


( 1 ) 



Initializes data that include thermodynamic constants in SI units 
LGNDRE; AS; ORDLEG; GAUSSG; GAUSSV 

calculate the Gaussian latitudes and weights and the associated legendre polynomial 
and tendencies using a rhomboidal truncation of wave number 15. 


ENTROROG; 

Interpolates 1° latitude x 1° longitude orography into 64 x 40 grid points that 
describes the finite difference horizontal resolution. 


INITIAL; ^ 

Reads the necessary data (sea-surface temperature and Orography, data sets) 
and prepares initial values to start the time integration. 

AMATRX; 

Calculates a matrix for solving K simultaneous equations for coefficients with given 
wave numbers. 


FAST; REALFT, FOUR1 
Perform fast Fourier transform. 


TR1; 

Transforms specific variable from physical space to spectral space. 

TR2; 

Transforms specific variable from spectral space to grid point space. 

RM1; 

Calculates the Fourier components using the spectral coefficients and the Associated 
legendre polynomials. 

RM2; 

Calculates meridional derivatives using spectral coefficients and Associated 
legendre polynomial tendencies. 

GUASS; 

Calculates the spectral coefficients. 

MMTRC; MTRCMX; 

Perform matrix multiplication. 

TR22; 

Calculates the meridional derivative in grid point space using the harmonic 
coefficients. 
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(3) . MODEL PHY SICS SUBPROGRAMS 

RADIATION; 

Calculates the longwave emissivities at the surface and the 5-level of model 
Atmospheres. 

REGR; 

Calculates second order regression coefficients of a series. 

INTRPL; 

Interpolates the zonal mean values at 5° latitudes into Gaussian latitudes. 

EMSIT; 

Uses the climatological 850-mb relative humidities to calculate emissivities at the 
surface and the 5-levels of the model atmosphere. 

SOLAR; 

Calculates surface and atmospheric solar radiation using a daily mean zenith angle. 
SURFTMP; 

Solves the thermal balance equation at the surface to calculate surface temperature 
and surface fluxes. 

CONVECT; 

Calculates the convective adjustment 
MADADJ; 

Performs the moist convective adjustment at a column of the 5-level model. 
CONVAD; 

Starts the convective adjustment scheme. 

CONDAD; 

Adjusts temperatures and mixing ratio of supersaturated levels with stable lapse 
rate below and above 

SMEQAD 

Called by MADADJ for the moist adjustment of temperature and mxing ratio 
MQADJ. 

Removes negative values of mixing ratio at a grid point using horizontal and vertical 
borrowing method. 

ESTABL; 

Calculates saturation vapor pressure 

(4) TIME INTEGRATION SUBPROGRAMS 
T22; 

Calculates explicit components of vertical temperature advection to apply the semi- 
implicit scheme 

RNUNV; 

Calculates the nonlinear dynamical terms 
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PREPARE; 

Calculates the hydrostatic geopotential at the five levels. 

VERTAD; 

Calculates the vertical advection 

SIGDOT; 

Calculates the vertical velocities 

TIME 1: 

Calculates the time-split vertical diffusion 

TIME 2: 

Calculates the semi-implicit time scheme 

TIME 3: 

Calculates the horizontal diffusion 

POTENG; 

Calculates the global potential energy and kinetic energy of the model’s atmosphere 
8.0 FY1991 WORK PLAN 

A 5-year simulation of the model's atmosphere will be established. The intraseasonal 
component will be isolated through the following operations: 

1 . Remove the climatological annual cycle by subtracting the 5-year mean value for each 
parameter. 

2. Transform the filtered data, for each year, to the Fourier space and remove the seasonal 
cycles (first and second harmonics). 

3 . Additionally (to the above low-pass filter) retain only amplitudes that have 30-60 day 
periods. 

4. Transform the truncated series to the physical space. 

Besides analyzing the differences in the intraseasonal variability between the different case studies 
(for example different sea-surface temperature patterns), we are interested in investigating the 
detailed relation between the thermal and dynamical patterns within each case. Excitation of the 
intraseasonal mode during El Nino episodes has been confirmed by many diagnostic studies using 
satellite data (e.g. Hanna and Crissman, 1991). The thermal heating in the model’s tropospheric 
column is a composite of the contributions of latent heat released during condensation of water 
vapor, radiational heating and the surface energy fluxes. At the same time, the tropical upper level 
divergence is, in part, a dynamical consequence to such thermal forcing. The model simulations 
will provide a continuous and comprehensive (in time and space) presentation of the thermal and 
dynamical fields. It is then possible to investigate, for example, how the upper level divergence 
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correlates with the tropospheric heating pattern within an intraseasonal time scale, and what are the 
spatial characteristics of the phase and amplitude of each parameter. 

The interaction of the intraseasonal oscillation with the local circulation and its effect on the 
medium and long-range weather forecast is not clear. To gain more insight into this problem, it is 
important to identify the development of the large-scale circulation associated with the low- 
frequency oscillation. The time filtering procedure described above will be used to extract the 
variability related to the low frequency oscillation of the temperature and geopotential height fields 
at the upper and lower levels of the model's atmosphere. The Empirical Orthogonal Function 
(EOF) analysis of the filtered fields will be applied for tropical and midlatitude zones. The eigen 
modes and eigen coefficients will be used to study the amplitude and phase of the time series. A 
reconstructed time series with the first few eigen modes can be compared with the monthly or 
seasonal statistics of a particular field. For example, Chen (1987) has used the EOF analyses to 
show that the first eigen mode of the filtered 200-mb temperatures is essential to form the warm 
center over the Tibet-India region and the cold center over the North Pacific during the Northern 
Hemisphere s umme r (1979). The low order eigen modes explain the largest variances of the time 
series. In this sense, if the reconstructed series using the first few eigen modes did not resemble 
the variance revealed by the infiltered time series that will lead to the conclusion that the variability 
at this time scale is not significant component of the large scale circulation. 

The modelling analyses will be linked and compared to analyses derived from satellite data 
such as the ISCCP-C1 data. Studies have identified the intraseasonal oscillation in the tropical 
OLR data, (Lau and Chan, 1985) and ER cloud temperature (Hanna and Crissman 1991). 
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Figure 2: A schematic representation of the sequence of operations for implementing one time step of the model. 
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